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Über erlegen; eines Schwinbers unter en Einfluß 
von Schwarz-Weiß-Regelungen. Fade 


ae "Herrn Prof. Dr. L. Prandtl zum 70. Geburtstag. 


_ Ergebnisse gemeinsamer Arbeiten von 1. Flügge-Lotz, H.F. ‚Hodapp, K. Klotter, H. Meisinger 
SE NAER Rs und K. Scholz. 


k re. / % Berichtet von Irmgard Flügge-Loiz, Berlin. 
f "Im der Technik werden aus Gründen apparativer Einfachheit oft unstetig Anbeitende ‚Regelungen 


so, daß bei ak Ausschlagen des Schwingers eine zusätzliche konstante Rückstellkraft oder ein kon- 
. stantes Rückstellmoment aufgebracht: wird, bei ‚negativen Ausschlagen dasselbe mit entgegengesetziem Vor- 
zeichen. Deshalb gehorcht die Bewegung in zwei aufeinanderfolgenden Intervallen verschiedenen, Gesetzen 
und muß stückweise berechnet werden, ähnlich wie Schwingungen mit Coulombscher Reibung. Das kostet 
viel Zeit, und es bestand deshalb ein Bedürfnis, Methoden zu entwickeln, dieeine Übersicht über den Be- 
wegungsablauf vermitteln, ohne diese stückweise Berechnung von Frequenz und Amplitude zu erfordern. 


In technies for the sake of simplieity of apparatus inconstantly working regulations are very often 
2; i constructed, in order to rectify swings of insufficvent self-dampening. These regulations f.i. are working 
Br z in such « way that with regard to positive turms .of the swing an additional constant force or, 
> momentum of reaction are produced; concerning negative turns the result is the same in converse order only. 
\ In consequence the motion of two successive intervals corresponds to different laws and has to be partially , 
calculated like the vibrations in conformity with Coulomb’s law of friction. T’his method, however, takes 
up much time; in consequence there has been a need to develop methods realizing a summary of the motional 

effect without requiring the partial calculations of frequency and amplitude. 


5 ‚ continue, destines ü ameliorer les syst&mes oscillants dont l’amortissement propre n'est pas suffisant. Ües 
- regulateurs travaillent p. exp. de telle fagon, que les £longations positives de l’oscillation produisent une 
force ou un moment reuctifs constants, tandisque les &longations negatives conduisent aux m&mes effects, 


mais de signe opposse. Par consequent, dans les deux intervall successifs le mouvement sera regi par deus . 
 lois differentes et doit ötre calcul& par fraclions, comme les oscillations avec le frottement de Coulomb. Cela , 


ewige beaucoup de temps et a conduit & developper les methodes nowvelles qui permeltent un apergu global. 
du mouvement sans avoir recours aux calculs fractionnes de frequence ou de l’amplitude. 


Aus yopomenna anmaparypbI B TeXHukKe yacTo IIPHMeHAIOTCA peryılatops IPepbI- 
BHCTOTO NeÜCTBUA, IEINBIO KOTOPEIX ABIACTCH YILYUIIeHNE KOMEÖATENBHEIX CHCTEM C HENOCTA- 
“ TOYHBIM COÖCTBEHHBIM 3aTyxaHuem. Ilonoönszie peryıaTopst Ppa60TaMT HAIPp. TAKHM O0Pasom, 
Fe; yTO HPH IIOJIOKHTEIBHOM OTKJIOHCHHH KONeÖATENBHOÄ CHCTEMEI HOABIAETCA NONOAIHHTENBHAA 
HOCTOAHHAA Peakıug MH PSAKTHBHEIH MOMEHT, & IPH OTPHNATEIBHOM OTKJIOHEHHH — TOIKE 
‘€ HPOTHBOHONO:KHEIM 3HAKOM. IloaToMmy, B AByX CieNyIIIMX APyT 34 APyToM UHTepBanax 
Korle0aHuAd HONINHAOTCH ABYMH PAsIHYHBIM 3aKOHAM U MOJIKHEIBEIYHCHATBCH IIO YACTAM, IIONOO- 
7 2.2 Ho KoleDaHusMm c Tpenuem KysoHa. ITo TpeöyerT MHOTO BpeMeHH H AAO IOBON K Pa3paborke 
HOBEHIX METONOB, IIO3BONAOILHX UPOCHeNHTB XON MBISKeHUA Öe3 HEOÖXONHMOCTH HONOÖHBIX 

APOÖHEIX BEIYHCHEHHÄ YACTOT HIIH AMINIHTYABI. 


h 1. Einleitung. In verschiedenen Zweigen der Technik treten bei Störung einer Hauptbewegung 
Schwingungen auf, deren Eigendämpfung den Forderungen nach schnellem Abklingen nicht 
genügt. Um in solchen Fällen rascher auf kleine Amplituden zu kommen, hat man. Regelungen 
geschaffen, die die Störbewegung möglichst schnell zum Verschwinden bringen. Diese Regelungen 
können stetig und unstetig arbeiten. Die letzteren werden als Schwarz-Weiß-Regelungen, auch 
als Ja-Nein-Regelungen bezeichnet. Über mehrere- bisher nicht veröffentlichte Arbeiten, die 
sich mit der Untersuchung von Schwarz-Weiß-Regelungen befassen, soll hier berichtet werden. 


2. Einläufiger Schwinger, Stellungszuordnung ohne Schaltverscehiebungen. Zunächst wurde 
der einläufige Schwinger untersucht, dessen freie (ungeregelte) Bewegung auf Grund der Gleichung 


RE RN Mr EA RE RE EIN TE BEER.) 
verläuft. Durch ein Stellglied greifen wir in den Krafthaushalt des Schwingers ein. Die Bewegungs- 
Be tung lautet dann 

+29 HW®o=Nß.. Yen (2), 
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gebaut, um Schwinger mit unzureichender Eigendämpfung zu verbessern. Diese Regelungen arbeiten z. B- 


Pour simphifier V’appareillage on emploit sowvent dans la technique des requlateurs 4 action. dis- 


BL Bß ist dabei 2 nahe des. ellgliedes # 
5 “ ‚stetigen Regelung wäre nun feine stetige Funktion von @ 0 | 
 Weiß-Regelung kann ß nur zwei ‚entgegengesetzt gleiche Werte 
Bereaeite. es NEHAuft: also nach den IIIGET 


Br: RR a N N a 
EN “ = Me 
Dabei sE zunächst ebene daß keine an re 


hr” 
RE durchgang der Regelfunktion P*—=x,9-+x;9 auch unmittelbar der Vorzeichenwechs 
N N ‘Deshalb Ener wir ‚statt ut Pr auch mit Es cr % 


| ea ER r 


NR e DE & * - Baht ex: = RR 
ee ch ER Be a N Re N FR 
A ‚arbeiten, wenn wir voraussetzen, daß. 8 ist. b = TE A 
Rt i Der Möglichkeit, dem positiven RE von nr ken eine 
re als auch eine positive Kraft zuzuordnen, entsprechen zwei 
‘verschiedene Regelsysteme (Bild 1). Ihren Sinn kann man in zwei . 

Grenzfällen leicht einsehen: 
“1. Wenn x sehr klein ist, wird us Regelfunktion rg, und 
N sen. F hat das gleiche Vorzeichen wie 9. Fassen wir nun das der zu- 
sätzlichen Kraft proportionale Glied NP, sgu F mit. dem ‚der ">. a 
ö Rückstellkraft proportionalen. Glied op zusammen, BI 
: 
’ 


at #+286 tar (p} ge =D Er da) 


+ Re: 


ER MERN aan / S0sehen wir, daß beim unteren Vorzeichen [d.h. dem er Pe 
pr und des! ern | Vorzeichen in Gl. (3)] im Mittel eine Verstärkung der Rückstellkraft 
bei Stellungszuordnung im System (System A), beim oberen (System B) aber eine Schwächung der Di 
eraud St Ne Rückstellkraft eintritt, die gegebenenfalls ‚den Schwinger labil 
' machen kann, 


2. Wenn x ; sehr groß ist, wird die Repeltunktiön F- x $,und sgn F hat ER gleiche Vorzeichen _ 
wie p. Fassen wir nun das der zusätzlichen Kraft proportionale Glied IEN Bo‘ sgnF mit dem der 
Reibungskraft proportionalen Glied 269 p zusammen 


#+20(p* Wassur) too EYE NILRR ae (4b), 


B so hen wir, daß das untere Vorzeichen (System A) im Mittel eine Verstärkung der Dämpfung 
Ba u.a hervorruft. 

JR Beide Grenzfälle zeigen, daß eine Regelung nach System A mechanisch günstiger sein wird 
als eine nach System B. Die genaue Rechnung wird diese Vermutung auch für andere x-Werte 
bestätigen. 

Die Lösung der Gl. (4) läßt sich zusammensetzen aus den Teillösungen 


a ea u Es ae 


dabei ist in jedem Intervall (Zeitraum, in dem F ein festes Vorzeichen hat) die Zeitzählung neu 

° zu beginnen und C, und: & sind die Konstanten, die im ersten Intervall die Anpassung an die 
gegebenen Anfangswerte gestatten. In den späteren Intervallen müssen die Werte C| und g, so 
bestimmt werden, daß ein stetiger Anschluß der Werte @ und 9 am Anfang des Intervalls an die 
Endwerte des vorhergehenden erfolgt. 


Die Funktion 9 (t) ist wegen des stetigen Anschlusses in den Schaltstellen in ihrem Gesamt- 
verlauf stetig. Da sich aber aus Gl. (4) ergibt, daß & (t) Sprünge an den Schaltstellen hat, muß $(t) 
‘dort Knicke haben. Daraus folgt, daß die Regelfunktion =p-+x ebenfalls Knicke an den 
Schaltstellen hat (vgl. Bild 1). 


Anstatt den Ausschlag 9 als Funktion der Zeit ‘anzugeben, kann man auch die Phasen- 
kurve (9) zeichnen, die einen guten Überblick über den Bewegungsverlauf gibt. Man kann nun 
zeigen, daß sich in einem passend gewählten schiefwinkligen 9—»lo- -Koordinatensystem 


<<. (9, plo) = arc cos = die Phasenkurve aus Stücken von kongruenten logarithmischen Spiralen 


on 1aßt die abwechselnd. Et Konvergenzpunkte ner = an =; ? haben. Der n 
s "ten Spiralenstücks hat inS ne na Ka 


Tun r - 6 
> BL 


EEE Se R De A IE ai 0) 0% ER Re) > 
' RE u 
£ nahe I rs ‚der Winkel zwischen Balkon r, und Se Radiusvektor am Anfang des. 
 Intervalls, der also vom Konvergenzpunkt zum Schaltpunkt Ss: = (9, 9|@) hinführt (vgl. Bild am 
ee! dieser Phasenebene. (9, lo) liegen alle Schaltpunkte F=9-+x09=0 auf einer Geraden, der 
Be „Schaltgeraden“, die die Phasenebene in zwei Hadbebenen teilt. Aus den einfachen Konstruktionen 


in der Phasenebene kann man die Eigenschaften der geregelten Bewegung leicht ermitteln. j Be 
’ 1. Im System A läuft in einem Intervall, d.h. in einem Zeitraum bestimmten Vorzeichens 
von F, die Spirale um den Konvergenzpunkt in der. URETDEIEN Halbebene, im Sal, DB um Te 
„.d dem. in derselben Halbebene (Bild 2 und 3): f N 
e r; y \ un! S‘ ; g ® A 
% n s; z i er ıfe | 


rare tg. xVw2=82 - 


T-xd4 
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Bild 2. Phasenkurve einer Bewegung, zusammengesetzt 
aus Stücken von eRSLuEzBi chen Spiralen; System 4. 


2. Im len Bist ein Verfehlen der Schaltgeraden möglich. Dann ruht die Regelung und h 
die Schwingung läuft wie eine ungeregelte (freie) um eine verschobene Mittellage (Bild 4). Man 4 
hat also durch die Regelung das Abklingen der Bewegung nicht verbessert. Um sich eine Übersicht | 
zu verschaffen, wann solches Ruhen bei den verschiedenen x-Werten auftritt, geht man am besten 
folgendermaßen vor: Da man den Bewegungsverlauf übersieht, wenn man die Folge der a 


} er % 


System B 


Bild 3. Phasenkurve einer Bewegung im System B; Bild 4. Phasenkurve einer Bewegung, bei der die 
7 gleiche Anfangswerte wie in Bild 2. Regelung zum Ruhen kommt. 
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punkte kennt; reicht es, wenn man auf jeder Schaltgeraden «=konst. die Schaltpunkte S, markiert. 
Man kann sich demnach einen Überblick über Bewegungsabläufe bei verschiedenem x verschaffen, 
wenn man in einer Ebene der Schaltwerte @,, 9,/®, die alle Schaltgeraden mit verschiedenem % 
enthält, Schaltpunkte einzeichnet. In dieser Ebene der Schaltwerte kann man nun Gebiete an- 
geben, die alle Schaltpunkte enthalten, nach denen Ruhen eintritt. | 


7 
System A 
9 System A 


TR 


Bild5. Endpunkt S. einer Bewegung. Bild 6. Phasenkurve einer periodischen Bewegung 
im System 4. 


3. Auf der Schaltgeraden F=0 gibt es eine Strecke mit besonderen Punkten. Wenn die 
Bewegung einen solchen Phasenpunkt erreicht, wird sie undefiniert, weil F zwar Null wird, aber 
kein Vorzeichenwechsel stattfinden kann (Bild 5 zeigt dies für System A). Dieses hängt mit dem 
geknickten Verlauf der Funktion F(t) an den Schaltstellen zusammen. Man kann für jeden Wert x 
diese Strecke ermitteln und dann ebenfalls in einer Ebene der Schaltwerte 9, und o,/® das Gebiet 
eingrenzen, in dem die Bewegung nicht mehr definiert ist. 
Solche Gebiete, Endpunktgebiete genannt, gibt es auch im 
System B. 


4. Es gibt periodische Lösungen, gegen die sich die Bewe- 
‘gung aufschaukelnd oder abklingend bewegt, wenn sie nicht auf 
Endpunkte oder Ruhepunkte zuläuft!) (Bild 6). 

In Wirklichkeit wird die Schaltung immer etwas nach dem 
Nulldurchgang der Regelfunktion F eintreten (Nachhinken des 
Stellgliedes), und deshalb tritt auch die Schwierigkeit der un- 
definierten Bewegung (Punkt 3) nicht auf. Nehmen wir an 
(vgl. Bild 5), die Schaltung erfolge immer um die Zeit ti, nach 
dem Nulldurchgang von F, dann folgt eine Bewegung hoher 
Frequenz (Bild 7), deren Mittelbewegung schnell auf 9—=0 und 


P=0 zuläuft, und zwar, wie Golling zuerst bemerkt hat, auf 
1 


‘ Grund der Gleichung g+x9=0 wiey=(;e * "Daraus 
folgt, daß man Regelungen mit kleinen positiven x-Werten bauen 
muß. Und da wegen der Möglichkeit des Ruhens der Regelung 
System B (vgl. Punkt 2) ausscheidet, ergibt sich für die Konstruk- 
tion einer brauchbaren Regelung die engere Forderung: Positiver 
Regelparameter x und System A. Weiter muß man bei einer 
Regelung wünschen, daß das Gebiet der Phasenebene, in dem 
die Schaltverschiebungen den Bewegungsverlauf wesentlich be- 


einf ist i i ; 
ee {lussen, groß ist in bezug auf die vorkommenden Werte PP 
kender Schaltung des. Stellgliedes weil man auf diese Weise leicht ein gutes Abklingen der Stör- 
im System 4; konstante Nach- bewegung erhält. Die Größe und Lage dieses Gebietes hängt von 


hinkzeit vi, = 5°. N Bo 
den Werten Re und x ab. Für verschiedene x-Werte übersieht 
man die Verhältnisse wieder am besten in der Ebene der Schaltwerte. Das Gebiet der Schalt- 


werte, nach denen die zugehörige Phasenkurve im Schaltstreifen verläuft, fällt dort im wesentlichen 
mit dem unter 2,3 definierten Endgebiet zusammen. 


ı) Für die mathematische Untersuchung dieser Frage vgl, Bilharz: Über eine gesteuerte eindimensionale Be- 
wegung. Z. angew, Math. Mech. 22 (1942), 206, 


ergibt u wa a ee en Stellungs- © 


MAN ua 
D 


rg a ae Rice Ne BER de Biken, ER 
Se) oder gegen einen Endpunkt und damit bei Schaltverschiebungen gegen eine hochfrequente, +1 Ma 
bei passend gewählten x(#> 0) im Mittel rasch auf null abklingende Bewegung; tu, 

| I oder wegen Eintritt des Ruhens der Regelung gegen eine a ar um an gegenüber 
Er Yan n Es der ungeregelten verschobenen Mittellage. 
RS; Einläufiger Schwinger, Stellungszuordnung mit a keerklebakene Die Wirkung von RR, 
_ Sehaltverschiebungen auf den einläufigen Schwinger mit demselben Regler wie früher wurde 
"genauer untersucht. Zur Vereinfachung ist die Dämpfung der ungeregelten Schwingung gleich 


2 N ? „N ull gesetzt. Das bedeutet für die Ergebnisse der Untersuchung keine wesentliche Einschränkung, 
9, »weil’ö'im. allgemeinen eine kleine Größe ist, denn sonst wäre keine Verbesserung der Dämpfung 
Be durch eine Regelung notwendig. ‚Für die Konstruktion bedeutet ö=0, daß die Phasenkurven- 


k stücke, die sonst im schiefwinkligen‘ | (9 »lo)-System logarithmische Spiralen sind, jetzt. u 
bogen im rechtwinkligen 9— plo-System sind, weil & (9; Plo) = arc cos 2 rj2 geht. 


0. Wir müssen zunächst den Begriff „‚Schaltverschiebungen“ erklären. Bislang erfolgebim 
" "Nulldurchgang von F* der Vorzeichenwechsel der Regelkraft, d.h. Kommando und Schaltung 
‚erfolgen gleichzeitig und zwar in dem Augenblick, in dem #* durch Null geht. Nun können 
sowohl Kommando als auch Schaltung unabhängig voneinander gegen den Null durchgang von 
F* verschoben sein. Folgende speziellen Fälle von Schaltverschiebungen wurden behandelt: 
- 1. Beim Nulldurchgang von F* wird ein Kommando gegeben, die Schaltung SR sobald || 
danach einen bestimmten Wert o, erreicht hat, und zwar erfolgt die 


Schaltung bei F*—= —o, für P*<0, - 
bei P*—= 1 0, für F*>0. 


4 x Bild 8 zeigt die Wirkung einer solchen Nachhinkzone mit konstantem 02. 
Ei 2. Beim Nulldurchgang von F* wird ein Kommando eben, ‚die Schaltung Sraelgr um die 
/ 


Zeit i, verzögert. 
system 4A Ri SystemB 


> 


| | F 
Bild 8. Verlauf der Regelfunktion F* und des Stellglied- Bild 9. ‚Verlauf der Regelfunktion #* und des Stellglied- ; 
ausschlags # bei konstanter Nachhinkzone o,. ausschlags $ bei konstanter Nachhinkzeit ,,. 


Bild 9 zeigt die Wirkung einer konstanten Nachhinkzeit, = 
3. Das Stellglied $ schlägt nicht von der einen Grenzlage + ßo sofort in die andere Grenz- 3 
lage + ß, um, sondern bleibt dazwischen solange- in der Nullage stehen, als |F*|<o, ist (tote : 
Zone). Kommando und Schaltung erfolgen gleichzeitig, sind aber gegen den Nulldurchgang 
: von F* verschoben. 


3 System A System B 
Ei KR Q: B=— Bo a P=+Bo 
- : EEE I Rune P=0 n>f*+ >20; P=0 
F*r< —0; ß=+ Po Pr< —0; B=—Bı 


Bild 10 zeigt eine tote Zone konstanter Breite o,: 


Bild 10. Verlauf der Regelfunktion 7* und des Stellgliedausschlags £# bei toter Zone o,. 


_  - vor, sondern können sich auch überlagern ; sie kommen. auch nich 

4; Form vor, aber die Klassifizierung erleichtert das Studium de 

Schlüsse zu ziehen über das Verhalten z. B. bei vorhandener, 

aM USW. LTR, Si N BEN ET ARN LEER e RE 
Tas Bakutstudiusn des Einflusses der Schaltverschiebungen müssen wii 

funktion #*=x,9 +x,9 benützen, weil die Schaltung nicht mehr beim DINIRIUPSHE 
erfolgt. Zur vollständigen Untersuchung müssen wir demnach zwei Parameter AM und E 


Die beschriebenen Typen von Schaltverschiebung = 


Dabei können wir wieder x, als positiv annehmen, denn durch das doppelte Vorzei e 
werden alle Möglichkeiten erfaßt. Da für den Fall, daß die Eigendämpfung Null is 
abschnittweise durch HIORCR, Rt IE h RS ni Be un 
| P—=20,c0s (+) H— mit D-aNBe 2... 


gegeben ist, hat die Steuerfunktion F* die Form 


P*=20, /»? + 0%} 005 (ot +e + u) + : 


Bild 11a. Phasenkurve Bild 11b. Ausschlag-Zeit-Diagramm 
einer Bewegung bei Regelung mit konstanter N achhinkzone im System A. 


Für eine konstante Nachhinkzone 0, ist in Bild 1la der Bewegungsverlauf für 
System A gezeigt. Da der Abstand von der Kommandogeraden +=0—=x,9-+x,9 dem 


Wert der Funktion F* proportional ist, ergibt sich, daß bei konstanter Nachhinkzone das 


Schalten auf Geraden erfolgt, die parallel zur Kommandogeraden sind. Wie bei der Re- 
gelung ohne Schaltverschiebung können wir auch hier leicht zeichnerisch die Bewegungsver- 
läufe ermitteln: 

1. Im System A verläuft für positive u die Bewegung abklingend.. Sie wird nach einer Anzahl 
von Intervallen, die von der Größe der Anfangsstörungen abhängt, in eine Bewegung, die nur 
zwischen den. Schaltgeraden verläuft, übergehen und schließlich (abgesehen von u=%0°) gegen 
eine periodische Bewegung um Pm= 90 streben. Die Schaltpunkte dieser periodischen Bewegung 


a 
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Sind die Schnittpunkte der Schaltgeraden mit der $/w-Achse (vgl. Bild 11a). Die Schwingungs- 
dauer 7 dieser Bewegung ist durch , 


j | Ws 


Bild 12a. Phasenkurve Bild 12b. Ausschlag-Zeit-Diagramm. 
einer Bewegung mit konstanter Nachhinkzone im System 4 bei # = 90°. (Die Regelfunktion ist der Geschwindigkeit 
proportional.) _ 


Bei u=90°, wenn die Steuerungsfunktion F der Funktion @ proportional ist, läuft die 
Bewegung ebenfalls auf eine periodische zu (vgl. Bild 12a und 12b). Diese periodische Bewegung 
setzt sich aber aus zwei ungleich langen Intervallen. zusammen. Die Schaltwerte 9, sind von 
Null verschieden. { 

2. Für negative u verläuft im System A die Bewegung angefacht (Bild 13). 


Bild 13. Phasenkurve einer angefachten Bewegung bei Regelung mit konstanter Nachhinkzone im System A. 


ei 
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3. Im System B verläuft die Bewegung für negative u gedämpft, bis Ruhen der Regelung i 
eintritt (Bild 14). Von da an setzt eine periodische Bewegung um einen im allgemeinen von Null 
verschiedenen Mittelwert ein. Wenn der ungesteuerte Schwinger eine Eigendämpfung hat (System- 
gleichung ö+269+@?p=0), tritt im Ruhegebiet eine gedämpfte freie Schwingung um eine- 


h Ber 


Regelung Regelung ruht 


arbeitet 


= - 


Bild 14. Phasenkurve bei konstanter Nachhinkzone Bild 15. Ausschlag-Zeit-Diagramm zur Phasenkurve 
im System B, bei der die Regelung zum Ruhen kommt. Bild 14. 


von Null verschiedene Mittellage auf. Die Gebiete, in denen Ruhen auftritt, können in einer 
Ebene der Schaltwerte (o,, 9,/®) für variierendes u angegeben werden. Bild 15 zeigt eine Bewegung, 
bei der Ruhen der Regelung auftritt, im Ausschlagzeitdiagramm. 

4. Für positives u verläuft im System B die Bewegung angefacht. 

Damit sind die Bewegungsverläufe bei konstanter Nachhinkzone untersucht, wenn wir u 
oder x,/x, variieren. Eine Variation von © = Yr?+ @?x2 bedeutet, daß sich der Abstand der Schalt- 


geraden = ändert. 


Tritt eine konstante Nachhinkzeit auf, so spaltet sich, wie bei der Nachhinkzone, die 
Schaltgerade in zwei Schaltgeraden g, und g, auf, die wieder symmetrisch zum Nullpunkt des 
Koordinatensystems liegen, nun aber im Gegensatz zu früher gegen die Kommandogerade um 
den Winkel ® t, geneigt sind (vgl. Bild 16). Dabei ergibt sich noch die merkwürdige Eigenschaft 
daß bei veränderlichem u, aber festen t, die Schaltgeraden g, alle durch einen festen Punkt 8 ; 
und die Geraden 9, alle durch einen festen Punkt S,, gehen. Die Koordinaten dieser Punkte ie 


b 
und $,, hängen nur von ot, und — ab. 
u} 


In Bild 16 ist eine Phasenkurve für Regelsystem A und u —o t, > 0 gezeichnet. Man sieht 
daß nach einigen Schaltungen die Phasenkurve nur noch innerhalb des Schaltstreifens verläuft; 
ein  dauerndes Weiterkonstruieren würde zeigen, daß sich.bei dieser Regelung (System en 
er bi N die en une en in nähert, deren Schaltpunkte die Schnittpunkte der 

chaltgeraden mit der p/o-Achse sind. Die Restschwingung hängt Sn . 
Schaltpunkte haben die Koordinaten SUNG MANBE Wessnhich Yon ans ele 


a b | cos (u —wt,) —cos u 
n o8 in (Lane ee (11) 
Die Schwingungsdauer 7’, ist durch 
vor [608 (u —ot,) —cos u 
1 Ta ss. .< (12) 


gegeben?). 


2) Für 4 — 90° ergibt si aus Tp =At 2 N i f 
RE [u 8 ch daraus 7’p 'y; Was zur experimentellen Bestimmung der Nachhinkzeit benutzt werden 


{ 


Bild 16. Phasenkurve einer Bewegung bei Regelung a 
Rate Na Sean ale im System 4; 4—0® bo >0. 


Y 


‚ Bild 17. Phasenkurve einer Bewegung wie in 
Bild 16, aber a <0. j 


H . 
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m Bild 17 ist eine geregelte ee im System A für ERnuLren u-Wert Kereitbuet, ‚sie 


. schaukelt sich auf, ganz entsprechend den Verhältnissen bei verschwindender Nachhinkzeit. Man 


‚kann sich leicht auch Beispiele für Bewegungen im System B bei verschiedenen u-Werten kon- 
'struieren. Da aber im System B wieder Ruhegebiete auftreten und die Bewegung je nach den 
Werten (u—ot,) entweder abklingend gegen die Ruhegebiete geht oder angefacht verläuft, ist 
das System B ungünstig. Für das bevorzugte System A kann man in der Ebene der Schaltwerte 
Grenzgebiete angeben, von deren Punkten aus die Bewegung in der Phasenebene immer nur im 
Schaltstreifen verläuft. Die Gebiete decken sich im wesentlichen?) mit den Gebieten, in denen 
bei der idealen Regelung ohne Schaltverschiebung die Bewegung undefiniert ist. Bei der 


Wahl des Regelparameters x =arc ig u, der wegen 


‚einer gut abklingenden Mittelbewegung klein sein soll, 


muß man darauf achten, daß wegen der Neigung der 
Schaltgeraden gegen die Kommandogerade der Win- 
kel (u—ot,) positiv bleiben muß, um im System A 
Abklingen zu erhalten. 

Die am Anfang dieses Abschnitts erwähnte Re- 
gelung mit einer toten Zone hat nicht die erfreulichen 
Eigenschaften der Regelung mit konstanter Nachhink- 
zone und -zeit. Die Gebiete undefinierter Bewegung 
verschwinden durch diese Art von Schaltverschiebung 
nicht. Bild 18 zeigt eine Phasenkurve für System A 
und positiven u-Wert. Sie setzt sich aus Kreisen um 


die Punkte p9= + 2 undp- O zusammen. Man kommt 


im System A bei u “, leichtzu Endpunkten. Das Grenz- 


gebiet in der Ebene der Schaltwerte, das die Endpunkte 
der Bewegung enthält, ist einfach zu bestimmen; wäh- 
rend sich bei der Bewegung ohne tote Zone Kreise 


(Eigendämpfung ö=0; wäre sie nicht 0, dann träten 


Ellipsen auf) ergeben, treten hier erweiterte Kardio- 
dien auf, allerdings nur für die Punkte a, und b,. Die 
gefährlichen Lagen der Punkte a, und b, sind etwas 
schwieriger zu beschreiben. Da die Konstruktion 


r 


Bild 18. Phasenkurve einer Bewegung mit 
toter Zone im System 4; u >. 


*) Geringe Unterschiede, weil Kommandogerade und Schaltgerade nicht wie dort zusammenfallen. 


\ 


dieser geregelten Bewegung aber Tele einfach. ist, 
nung sich einen vollen. Jberblick über den Ablauf ver 
. dämpfung d=0 eine tote Zone unbedingt zu vermeiden, da sie 
gegenüber der idealen Regelung bringen kann und gar keine Vor Zi 
mit sich bringt. Ist dies nicht ‚möglich, ‚so ‚muß man dafür Sorge trage 
_ hinken der Schaltung die Bewegung immer definiert bleibt. und si 
periodische Bewegung nahe am Nullpunkt der Phasenebene geht, 
0 Wenn man den Einfluß von Schaltverschiebungen ‚auf die geregelte wegung von 
gern mit Eigendämpfung ö untersucht, Sun: Big! geringe Modifikationen en Resulta 
ö=0, die leicht zu ermitteln sind. SR TE 


Zur! 


> Er 

IE 4. Einläufiger Schwinger, Hattessch indie ohne ei: mit Rückführung, Die: : 
bisher behandelte Regelung arbeitet mit Stellungszuordnung, d.h. der Ausschlag des Stellgliedes 
hängt vom Ausschlag, p und der Geschwindigkeit 9 ab. Vielfach. werden nun auch Rege- “ 
lungen gebaut, bei denen die Geschwindigkeit ß des Stellgliedes von den Werten pund abhängt. re 
Solche Regelungen bezeichnet man als Re gelungen mit Laufgeschwindigkeitszuordnung. 
Wenn die Geschwindigkeit des Stellgliedes nur zwei feste entgegengesetzt gleiche Werte 2 

‘ annehmen kann, sprechen 'wir von einer Schwarz-Weiß-Regelung mit Laufgeschwindigkeitszu- Fr 
ordnung. Sie unterscheidet sich von der Schwarz-Weiß- -Regelung mit Stellungszuordnung 
wesentlich dadurch, daß das eo nicht nur zwei REN einnehmen Ar en BE 
‚seine Lage stetig ändert. / Er 


- Für einen einläufigen 'Schwinger eds die. geregelte Bewegung bei Laufgeschwindigkeits- NY ee 
. zwordnung mit und ohne Rückführung untersucht, we keine Rückführung vorhanden ist, 
Ki lauten au System- und Beaelelchunn SASSFILSE | ae 

Be +29 +0 NB.. ne ee 1 2 BE Var 


Erz = + Vogur- zur Mr pt rap IR a en . (14). Seo): = 


Ab t 


Die beiden Vorzlichen ne eheider zwei Bean die in Analogie zur Festsetzung i in. 2 
Abschnitt 2 als A und B bezeichnet sind: - 


B= —VsgnF+ system A. 1% (15) 
ß=+VsgnF* _ System B : 


Man kann hier aber nicht den Grenz- 
fällen = = x —0 und x—» eine ein- 
H 


1 3 
fache mechanische Deutung geben. In 
Bild 19 ist für System A und positives x 
ein angefachter Bewegungsverlauf gezeich- 
net. Eine Regelung mit Stellungszuordnung- 
hätte in diesemFall (4,=2>0) zu einer ge- 
dämpften Bewegung geführt. Ersetzt man 
den in Bild 19 angegebenen ß-Verlauf durch 
einen stückweise konstanten im Sinne einer 
Stellungszuordnung (punktiert eingetragen), 
so erkennt man, daß #< 0 sein muß, wenn 
=P-x9 anden neuen Schaltstellen ver- 
schwinden soll.‘ Dem entspricht aber bei 
Stellungszuordnung F eine angefachte Be- 
wegung. Allerdings bestehen nicht immer 
so einfache Analogien zwischen Regelun- 
gen mit Laufgeschwindigkeitszuordnung _ 
ohne Rückführung und  Stellungszuord- 
nung. Eine wesentliche Rolle spielt näm- 
lich der Anfangswert des Stellgliedaus- 
schlags ß.' Dieser ist in dem Bild 19 zu- 
grunde liegenden Beispiel Null. Ist er 
h a a et wesentlich von Null verschieden, wie etwa 
matt Laufgeschwinäigkeitäterelung ulne Hhekiimene an Bild2i, dann ee 
System 4; Regelparameter @ x = 0,1, Eigendämpfung D=0. der Ausschlag des Stellgliedes während des 


Pl /w 


einer ee Die we einer 
Ds Aalen: Kine eine ig führt also nicht etwa nur 
er ‚die Bewegungsvorgänge BULES. zu ‚machen, sondern. ‚kann auch Ben andere Be- 
. wegun, stypen ergeben.“ , De 


Br { # x ER ie analytische Lösung für die Teilintervalle ist leicht anzuschreiben, v wenn man beachtet, 
ey! daß die einmal differenzierte di. 13) DEE REF | E Ri 

} AR Re: Ne: VE, \ f: Bi I I > ri ‘ 
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- Bild 20. Ausschlag-Zeit-Diagramm einer Bewegung mit Laufgeschwindigkeitsregelung ohne Rückführung . 
im a, veleny 4; ox=5b,D=0,1. Da der zugehörige periodische Schaltpunkt im Koordinatenursprung liegt, 
p ; klingt die Beprogung nach 0 ab. 


©, eund Bo sind durch die Werte von @, p und ß am Intervallanfang bestimmt. In Bild 20 ist der 
Verlauf einer Bewegung im System A bei positivem u-Wert durch die Diagramme g(t), &(t), Fi), 
40) und ß(t) dargestellt. Aus der Diskussion der Gl. (16) ergibt sich, daß'g (t) Sprünge hat, also 
o(t) Knicke. Deshalb hat im Gegensatz zur Regelung mit Stellungszuordnung die Regelfunktion 
eine stetige Ableitung F*—=x,9-+x,9. Daraus folgt, daß bei der Regelung mit Laufgeschwindig- 
keitszuordnung keine Schaltpunkte auftreten, in denen die Bewegung undefiniert wird. 

Die Phasenkurve der Bewegung ist durch die drei Koordinaten o(tf), p(t) und £(t) t) bestimmt, 
ist also eine Raumkurve. Sie setzt sich aus Stücken ‚‚auseinandergezogener Spiralen‘‘ zusammen. 
Eine so einfache Konstruktion wie im Fall der Stellungszuordnung ist nicht mehr möglich. Für 

- "einen Überblick aber den Verlauf der Bewegungen reicht es, wenn man die Folge der Schaltpunkte 


Bild21a Anfangs- und Endisochronen 
Y für eıne mit Loufgeschwindigkeitszuordnung ohne 
Rückführung gesteuerte Bewegung. . 
System A wx=1 D=-41 
- Änfangsisochrenen (Intersellange z 8 30) 
N — Endisochronen (Intervallonge 2.8.\3 ). 
Die Bezıfferung anden Punkten —e— der Endisochronen weist hin af die zugehörigen Schallaerte g, der Anfangsi: 3 


Bild 21a. Anfangs- und Endisochronen für eine mit Laufgeschwindigkeitszuordnung ohne Rückführung geregelte Bewegung. 
System A; ve =1, D = 0,1; periodischer Punkt P. N ® A 


Anfangsisochronen genannt. Ebenso kann man nach vi,cna bezifferte Endisochronen zeichnen, 
die alle Schaltpunkte 5 zusammenfassen, in denen Intervalle gleicher Länge »tjend enden. 
Längs jeder Endisochrone läßt sich noch eine Bezifferung anbringen, die angibt, mit wel- 
chem Schaltwert , das hier endende Intervall angefangen hat. Ergänzt man dieses Bild 
noch durch die Anfangsisochronen der Schaltpunkte b und die Endisochronen der Schalt- 

’ punkte a, so hat man eine Kurven- 
tafel, die es erlaubt, zu einem be- 
kannten Schaltpunkt s„(9,,ß,) die 
Bestimmungsstücke des darauf fol- 
genden %&S,+, zu bestimmen. In 
Bild 2la ist eine solche. Schalt- 
punktiolge eingetragen. Die Be- 
wegung beginnt mit dem  Schalt- 


N 
punkt So — (or ß; 2) anna (8.9, —5, 8) 


im 2. Quadranten. Er liegt auf der 
Anfangsisochrone t,=210° und ist 
ein Schaltpunkt vom Typ a. Der fol- 
| ende Schaltpunkt 8, liegt also auf 
Bild 21b. Ausschlag-Zeit-Diagramm zu der Bewegung, deren der Endisochrone vom Typ 5b mit 

Schaltpunkte S,, in Bild 21a eingezeichnet sind. Vijena—210°. Diese läuft durch den 


*) Für die Definition der Schaltpunkte a und b vgl. Bildi1. Im Schaltpunkt a springt bei Stellungszuordnung der 
Ausschlag des Stellgliedes, bei Laufgeschwindigkeitszuordnung die Ges chwindigkeit des Stellgliedes von seinem negativen 
Wert auf seinen positiven, im Schaltpunkt b umgekehrt. } E 


Be 5 


den adraten. ee anne leicht zu nd Punkten‘ Faser, ie nn N 
müssen wir feststellen, auf welcher Anfangsisochrone vom Typ d der Punkt 8, liegt. Durch Inter- * 
- polation zwischen den Anfangsisochronen db mit 9 tan = 120° und 140° finden wir » tyant 138 
‚ und suchen nun den Punkt S,. Er muß auf der Endisochrone vom Typ a mit der Ziffer 1310 


diesem Fall 


. neue Steuerfunktion Gl. en bestimmt ist. 
(Bild 22). 2 


Stellen. Deshalb hat auch die Regelfunk- 


‚Ursache für das Auftreten von Schalt- 


liegen usw. In Bild en ist ‚das zu der eben a Se Ban Ausschlagzeit- 


 diagramm gezeichnet. en N 


Im System A ist Er Kmebeschwiniigken‘ Er gering. Liegt jedoch ein periodischer 


m | Punkt in der Schaltebene, so strebt die Bewegung diesem zu. Dabei kann trotz eines kleinen Be 
> Schaltwertes @, der Maximalausschlag von p noch groß sein. Im System B salat die Bewegung. 
nach wenigen Schaltungen bei Ruhepunkten. 


Nach diesen Ergebnissen empfiehlt sich eine ER mit en 
ohne Rückführung nicht. Ganz anders liegen dagegen die Verhältnisse bei einer Regelung mit 
Laufgeschwindigkeitszuordnung mit ee ‚Die Regelfunktion Er lautet in 


RG | munptmptmäß. AT AT N .. (2l). 


Die Tösinsen lassen sich für die einzelnen Intervalle in genau derselben Form wie bei der ee 


ohne Rückführung Gl. (17, 18, 20) anschreiben. Natürlich ist auch bei gleichen Anfangswerten 


wie früher (Bild 20) der Verlauf der Bewegung anders, weil die Länge der Intervalle durch die 


Da die Geschwindigkeit: ß des Stell- 
gliedes an den Schaltstellen unstetig ist, 
hat der Ausschlag ß() Knicke an diesen 


tion F* ()=rXı 9 + p-+%ß N wieder 


9% 
Knicke an den Schaltstellen wie bei der 
Regelung mit Stellungszuordnung. Da, wie 
wir früher schon bemerkten, Knicke der 
Regelfunktion F* an den Schaltstellen die 


punkten ist, in denen die Bewegung un- 
definiert wird, werden im Raum der 
Schaltwerte o,, 9,/o, ß nun wieder Gebiete 
auftreten, die solche Punkte enthalten 
(Randpunktgebiete, kurz Grenzgebiete ge- 
nannt). N 
Bei beliebig vorgegebenen Anfangs- 
werten kann man die Bewegung am einfach- 
sten wieder wie früher durch die Folge der ; Bild 22. Ausschlag-Zeit-Diagramm einer Bewegung bei 
Schaltpunkte in der Schaltebene mit ein- Laufgeschwindigkeitsregelung mit Rückführung im Sy- 
gezeichneten Anfangs- und Endisochronen ee a N 
verfolgen. Die Mannigfaltigkeit der Schalt- RESET 

ebenen ist größer als bei der Regelung 

ohne Rückführung, weil zwei Parameter x,/x,; und x3/x, die Lage der Ebene bestimmen. Die 
Schaltebene enthält den Nullpunkt des Phasenraumes, liegt aber im Gegensatz zu früher schräg 
gegen das Achsenkreuz, so daß sie keine der Achsen enthält. Bild 23 zeigt eine Projektion einer 
Schaltebene mit Anfangs- und Endisochronen in die @,,ß,-Ebene. In der Projektion der Schalt- 
ebene wird auch das Gebiet der Schaltpunkte sichtbar, in denen die Bewegung, wenn keine Schalt- 
verschiebungen vorhanden sind, undefiniert wird (Parallelstreifen, in dem nur Endisochronen 
vorhanden sind). Da in Wirklichkeit aber Schaltverschiebungen immer vorhanden sind, wird 
man wieder eine hochfrequente Bewegung erhalten, deren Mittelbewegung leicht bestimmt 
Er kann. Für sie gelten die Gleichungen 


ee ee (22) 


> 


von G1.(23) in Gl.(22) ergibt sich . 


Äh Durch Einsetzen 


Es 


TRITT N 


irrt... 
Ra EC 
— Das heißt aber, daß die 1 
wieder ‚eine freie Schwingung ist 
Dämpfung und Rückstellkraft ab 
über der ungeregelten Bewegung 
sind. Die Schwarz-Weiß-Regelung 
also nach Erreichen eines Schaltpunktes 
im Grenzgebiet infolge der Schaltverschie- 
bungen wie eine stetige, Regelung. Ds 
ist eine wichtige Erkenntnis. Sie bedeutet, Bo 
‘daß man mit der apparativ primitivrn 
Schwarz-Weiß-Regelung. ähnliche Erfolge 
wie mit einer stetigen Regelung haben 
kann, vorausgesetzt, daß man mit geeig- 
neten Schaltverschiebungen im Grenzgebiet 
arbeitet. Da die Regelung nur Sinn hat, 
wenn die geregelte Bewegung besser b- 
klingt als die ‚ungeregelte, ergeben ich 
aus Gl. (24) Vorschriften über die Wahl dr y 


 Regelparameter und Ss Bild 24 zeigt das 
1. 1 
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er Dämpfungsmaß | 


63 Tre x 
ayı+? Vı+ wie 
. Kg } ; Ka pre 


: 2 en | ‘wenn D=Öö/® das Dämpfungsmaß der un- 
| System A:wx,=23, 2514 | geregelten Bewegung bedeutet. 
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Bild 23. Anfangs- und Endisochronen für eine Bewegung bei _ f ER 
Laufgeschwindigkeitsregelung mit Rückführungim System 4; - 5. Dreiläufiger Schwinger, Stellungs- 


Be AT zuordnung. Nachdem die Eigenschaften 

des geregelten einläufigen Schwingers aus- 

führlich untersucht sind, wenden wir uns nun der Regelung eines mehrläufigen Schwingers zu, 
z.B. der eines Flugzeugs. x 

4 Ein dynamisch stabiles Flugzeug durchfliege eine gradlinige Bahn, die durch den Anstell- 

h 7 \ winkel &,, den Längsneigungswinkel 9, und die Flugzeuggeschwindigkeit ®, gekennzeichnet ist. 
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Bild 24. Dämpfungsmaß D der Mittelbewegung bei Laufgeschwindigkeitsregelung mit Rückführung. 


A ü N Sobele es in seiner Bahn gestört wird, führt e es eine D eedamarte Beie Schwingung aus. Die alre da 
dieser Schwingung auftretenden Abweichungen des Anstellwinkels, des Längsneigungswinkels und 


RN die i im allgemeinen geringe Dämpfung der Bahnschwingung zu vergrößern, sei eine Regelung ein- 
N gebaut. Die Lagenabweichung. ® und die Drehgeschwindigkeit Ö werden durch Kreisel gefühlt, 


und dadurch das Höhenruder so gesteuert, daß die entsprechende Auftriebs- oder Momenten- : 


änderung die Bewegung beeinflußt, d.h. möglichst zusätzlich dämpft. 3 


\ 


AR BE Wenniwir. z.B, das‘, Höhenruder regeln und annehmen, daß das Gleichgewicht in, > 


E72” Auftriebs- und Eahggiehkune} Posch nicht beeinflußt wird, lauten ‚die Sem und Regel. 
E Belastungen 5 


t & . + . { ’ nt 


» | d Es N re 
anetbutt lite, ze et Ba A, 


en e = a Re “ ne (es) 
. em). 


nn tms 0 +0. Re RA 


T bezeichnet die dimensionslos gemachte Zeit t. 


Da ng während eines Intervalles konstant ist und die Systemgleichungen lineare Differential- 


die Form: 


I= EB, Pi 7 Bd; = IA, ge au Ur 30h gucde 4. C. . (08). 

Er Dabei DS ichnen die Konstanten 4, B, € die durch die Regelung hervorgerufenen Mittellagen: 

verschiebungen, die ®; werden durch die Anfangsbedingungen bestimmt, die C; und A; sind durch 
die Konstanten der Systemgleichungen bestimmte Proportionalitätszahlen, die A; sind die Wurzeln 
der charakteristischen Gleichung des Systems. Sie hat im allgemeinen zwei Paare konjugiert 

(] 

3 komplexer Wurzeln. Demnach können wir vier Anfangswerte y,— — In du da ER und 13 
erfüllen®). Genau wie früher wird sich zeigen, daß der Verlauf der geregelten Bewegung wesentlich 
vom gewählten Regelsystem und von den Anfangsstörungen abhängt. Die Darstellung macht 
allerdings einige Mühe, weil man wegen der vier Anfangswerte eigentlich im vierdimensionalen 
Raum arbeiten müßte. Wenn man sich aber auf die Darstellung der Schaltwerte für festes x be- 
schränkt, kann man in einem dreidimensionalen Raum arbeiten, ‚weil 9 und ö' dann Se mitein- 
ander durch die Beziehung 


tr. (29) 


V9 + 9, ? gewählt. 


Nach den Ergebnissen der Abschnitte 2 und 3 wird man versuchen, die Regelung möglichst 
so zu bauen, daß man nach wenigen Regeltakten auf Schaltpunkte im Grenzgebiet stößt, nach denen 
die Störung in schnellen Schwingungen rasch abklingt. Deshalb hat es einen Sinn, sofort die Frage 
nach diesem Gebiet anzuschneiden. 


verbunden sind. Als Koordinaten dieses Schaltraums®) sind y v und &= 


- \ ; 
: 5) Die Systemgleichungen schreiben sich am bequemsten in &, 9 und v. Die Störungen sind aber meist in y,. 9, 9° 
und v gegeben, wobei y = d — x den Bahnwinkel bezeichnet. 
°) Zur Gedächtnisstütze seien die Bezeichnungen noch einmal zusammengestellt: 
1. Einläufiger Schwinger mit Stellungszuordnung 


Ebene der Schaltwerte: Koordinaten 9, und Pg Schaltstrahl oder Schaltgerade: 9, = —% P2. 
2. Einläufiger Schwinger mit Laufgeschwindigkeitszuordnung Ä 
Raum der Schaltwerte: Koordinaten 9, 9; und ß. Schaltebene: Koordinaten 9, = —* go, und ß. 
3. Dreiläufiger Schwinger mit Stellungszuordnung x 
Raum der Schaltwerte: Koordinaten v5, Yyg 9 ds’ (4 Dimensionen). 
4 Schaltraum: Koordinaten v,, Y3 9 = —xg ds’ (3 Dimensionen). 2 


Grenzgebiete (oder Endpunktgebiete bei fehlenden Schaltverschiebungen) sind in der Ebene bzw. den Räumen der Schalt- 
werte. eingezeichnet, sie werden in den Gebilden geringerer Ordnung (Schaltstrahl, Schaltebene und Schaltraum) im Schnitt 
sichtbar. - 


der: Geschwindigkeit von den stationären Werten sollen mit.«, ® und » bezeichnet werden. Um = 


gleichungen mit konstanten Koeffizienten sind, haben ZRISSheH zwei Schaltungen die Er 
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112 / Flügge-Lo N 2, Über Bereguisen eines Schwingers 


Zunächst stellen wir fest, daß die Regelfunktion wieder Knicke aufweist, die, wie früher 
angegeben, für das Auftreten von Grenzgebieten im Raum der Schaltwerte verantwortlich sind. 
Aus Gl. (26, III) folgt, daß 9’ Sprünge hat, also 9 Knicke, damit hat auch F=d9+x%' Knicke. 
Die Berandung des Grenzgebiets ist durch die gleichzeitig geltenden Bedingungen #=0 und F=0 
gegeben. Im oben definierten Schaltraum erfüllt das Grenzgebiet den Raum zwischen zwei parallelen 
Ebenen. Die Größe des Grenzgebietes hängt vom Regelparameter x und von den Beiwerten der 
Systemgleichungen ab. _ Be Be) 

* Die Mittelbewegung nach dem ersten Schalttakt im Grenzgebiet ist durch die Systemglei- 
chungen (26) und die Gleichung IS 


BAM TI Te er Bee (30) 
bestimmt. Um die Lösungen o,,, Os Yms HR zu finden, werden wir die charakteristische Gleichung 
u) di Co+ E14 0 
Mt d bunt bh: C20 Dr Me . (31) 
At AA dit 634° 0° N 
0 1+x4 0 0 
auflösen müssen. Man sieht sofort, daß sich die Determinante in Gl. (31) reduziert; 
i a Co + ey A 
(1+»4 10 10 11 Bein“ 
Ag Ayıd 27) 


Es ergeben sich drei Wurzeln 
1 
em, } 
1 2 /, und A, 


A, und A, hängen nur von den aerodynamischen Beiwerten ab und sind bei einer Regelung über das 
Höhenruder negativ reell. Die Wurzel A| wird negativ, wenn x positiv gewählt wird. 


Um die Wirkung einer Schwarz-Weiß-Regelung zu zeigen, fügen wir zwei Bilder an (Bild 25 
und 26). In Bild 25 wird die Wirkung einer Regelung mit Stellungszuordnung auf die Variabeln 9, d° 


Bild 25. Dreiläufiger Schwinger, geregelt mit Stellungszuordnung im System A 
Regelparameter # =5,. Die Bewegung führt ins Grenzgebiet und ist dort als 
Mittelbewegung gezeichnet. 


und v/v, gezeigt, wenn die Schwingung durch eine Anfangsstörung der Geschwindigkeit hervor- 

gerufen ist. Im Bild 26 ist die zugehörige Bahn-Abweichung aufgezeichnet. Zum Vergleich sind 

die geregelte Bewegung mit kleinerem x-Wert und die freie Schwingung eingetragen. Entsprechend 

den Ergebnissen des 3. Abschnittes ist die nach dem Endpunkt eingezeichnete Bahnabweichung 

die „mittlere“ Bahnabweichung. Die wirkliche Bahn ist eine schnelle Schwingung kleiner Amplitude 

u Bun Ihre Frequenz und Amplitude hängen von der Größe der Schaltverschie- 
ungen ab. 


Bei Be Kar man für den en See ebenalls, die 3 
ee der Grenzgebiete bei Regelung mit Rückführung durchführen. In diesen wird dann 
wieder eine Mittelbewegung aa u einer stetigen Steuerung entspricht. 
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Bild 26. Abweichungen von der stationären Flugbahn beim Beispiel des Schwingers von Bild 25. Vergleich mit u 
Bahnabweichungen bei Regelung, mitx = 0,1 en bei ungeregelter Bewegung unter den gleichen Anfangsbedingungen. 


Del, 
SZ = er Zusammenfassung. 

BR, DIE) Untersuchung über Schwarz-Weiß- „Regelung haben: gezeigt, daß es rer ist, 
die Regelung so zu bauen, daß man nach wenigen Regeltakten auf Schaltpunkte im Grenzgebiet 
(Randpunktgebiet) stößt, nach denen die Bewegung zu einer Kriechbewegung wird mit.überlagerter 
= - schneller Schwingung kleiner Amplitude. Die Ausdehnung dieses Grenzgebiets in der Ebene oder 
dem Raum der Schaltwerte (je.nach der Ordnung des Problems) hängt von den Koeffizienten der 
Systemgleichungen und den Regelparametern ab. Das Abklingen der Kriechbewegung hängt 
ebenfalls von den Regelparametern und von Beiwerten der Systemgleichungen ab, die Rest- 
schwingung um die Kriechbewegung wird durch die Größe der Schaltverschiebungen bestimmt. 
Diese müssen deshalb klein gehalten werden. 


a a a a Zn Ba = ae un u u 
Normulobweichung 


Theorie und Praxis der Periodogrammanalyse 
von Registrierkurven, die im wesentlichen aus | er 
nichtpersistenten Wellenzügen bestehen. a 


Er. : a | | Von H. Blume. 


Die mathematischen Bedingungen der Quasipersistenz werden für den Fall abgeleitet, daß die Registrier- Ve 
kurve im wesentlichen aus Wellenzügen konstanter Frequenz und Phase, aber beliebig veränderlicher Am- a 
Be besteht. 5 
"The mathematic conditions of the ““Quasi- „persistence” (= persistence so to speak) are demonstrated 
PR the case supposing that the registering curve is essentially consisting of undulations of constant frequency 
and phase, but of arbitrarily variable amplitude. 
+ “ Des conditions mathematiques de la quasi-persistance sont demonirees pour le cas ou la courbe en- 
S registree_est essentiellement constituee par un lrain d’oscilations de phase et de frequence constantes, i 
mais d’une amplitude variable de fagon arbitraire. 


Maremarnyeckue yYCJHIOBUA ‚‚CBEPX-YCTOHYUBOCTH‘‘ (KBASHNEP3HCTEHN) BEIBONATCH MIA 


Toro cHyyaA, KOTNA PEeTHCTPAIHOHHAF- EPHBAA COCTOHT u3 pAna Kolebannf . NOCTOAHHOH 
yYaCTOTEHI uU da3bl, HO IIPOHSBOJIBHO NSMEHAIOIMEÄCH AMIIIUTYAH. 


Die praktische Verwendbarkeit der Periodogrammanalyse beruht auf ihrer Tragweite für 
die Frequenz- und Amplitudenberechnung über den einfachen Fall einer persistenten Welle - F 
hinaus, aus dem sie abgeleitet ist *).. Denn in der Praxis kommt es kaum vor, daß eine Registrier- 
kurve wesentlich aus persistenten Wellenzügen besteht. Es ist deshalb grundsätzlich wichtig, 


1) Stumpff: Grundlagen und Methoden der Periodenforschung, Berlin 1937; Blumeund Kolb: Über die 
Analyse von Registrierkurven von Schwingungsmessern mit Hilfe des Harmonischen Analysators nach dem Periodogramm- 
verfahren. Z. Instrumentenkde. (noch nicht erschienen). 
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s 1. Periodogtammanalyse ei einer Kun "konstanter Frequenz ı 
“z U veränderlicher Amplitude. REN 
neh sei ein Wellenzug der analytischen Form y= (at+b). cos a 


Das Analysenintervall der Länge p soll sich vom Punkte g> 0 bis q+ p.der waagerechten “a 
ö) ‚erstrecken, ı und die ra der Kurve soll so es werden, daß De get er ER 
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‚Bild 1. Kurve kon Gen Era mit Eines Bereich konstanter und linear ae Se > ie. 
Ä £ abfallender Maapgendett: ; De E 5 s 
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weise um Re ten Größe A x 0 verändert wird. Für einen ee Pünkg q Yet. SR Et 
dann der Harmonische Analysator die Lösung des verallgemeinerten Fourierkoeffizienten, dern 


in komplexer Schreibweise lautet: a | oe Br 
qa+P» ALER N 1 ar» 935 ; PAR. % 
a(x, g+ib(e, = [ur b)eils-mt—astAds+ N (a4 5) er 2. u 
, q q 4 se eh 
Die Auswertung: ergibt nach einiger Umformung: = = Kr: 2 
a(a,g)t+ib(2,g)=(ag+b) = - PR ine p2—ag+ M 
— a) pj2 | 
„finale Neo —an+ m —eite—nn nat) £ a 
x —.a (2 —o) p]2 er 
sin («+ a) p/2 RrtIoN 
a Di en eilte +) pl + a0 —M 
+ (ag+b) en | 
RAR an Re PR eiliat apa +aa— Bl — eilt Stel 
en (x +) pl2 


=ys,)+n(a,gq 


A Die, Summe der beiden ersten Integrale in Gl. (1a) bildet den Häuptyekter. Sie sei mit 
BEN, y(x, q) bezeichnet, die Summe der beiden letzten der Störvektor 7 (x, g). 

Be. Es werde zunächst gefragt, welchen Wert der Hauptvektor für den Fall = « annimmt. 
Die Durchführung des Grenzübergangs liefert rechnerisch den Wert: 


? 110 0= (an +54 Lap)eimsetn Be RR ee 


\ Das Resultät bedeutet folgendes: Stimmt die Versuchsfrequenz &, = Een (n=11/2..) für 
i Haie, 2 e } n ? 
ein a=n' mit der Frequenz « = RT, überein, (n =r) so besitzt das Phasendiagramm für 


diese Versuchsfrequenz x, die Neigung —&. Die Amplitude in der runden Klammer stellt das 
arithmetische Mittel der Ordinaten der linearen Gleichung 2= at + b innerhalb des Analysen- 
intervalles » dar, wie Bild 1 veranschaulicht. Infolgedessen zeigt das Amplituden-Diagramm 


’)H. Blume: Über die Analyse kurzer Kurvenzüge. Z. angew. Math. Mech. Bd. 23 (1943), S. 346358. 
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| tie und Praxis der Periodogrammana 


A Be rar 3 5 LT: ne DEZ Er RE PEN er han: 
y(o, g) in Abhängigkeit von q die Veränderung dieses Mittelwertes, d.h. es ist e 
# ‚derselben Neigung, wie der lineare Abfall der Amplitude der Ausgangswell.‘ * 
Man kann die Bedingungen feststellen, unter denen der Störvektor n(x, q) für den Fall 
%© = & gegenüber dem Hauptvektor vernachlässigbar klein ist. Die Abschätzung lehrt, daß der 


Absolutbetrag Hedi 2- 7; wenn man ein Glied — 127 vernachlässigt. Dabei. 


r(n-+-r) 


m(n+ r)* 


gibt, die Zahlr an, wie oft die Wellenlänge der gegebenen Frequenz « in dem Analysenintervall 


enthalten ist, wobei r eine positive rationale Zahlist. Für r = 6 ist dieser Betrag kleiner als !/3.. 


Setzt man diese Annahme als erfüllt voraus, so werden also das Amplituden- und Phasendia- 
gramm von Gl. (1a) für den Fall % = a im wesentlichen durch den Hauptvektor allein 


dargestellt. 


Es bleibt noch zu zeigen, daß der Absolutbetrag des Hauptvektors y(e, g) nur in der 


Umgebung von &—=& wesentlich von Null verschieden ist. Von dem ersten Summanden in. i 
‚Gl. (1a) ist es bekannt, daß er sein größtes Maximum mit dem Wert |ag+blin @= «besitzt. 


Es bleibt also zu beweisen, daß dies auch von dem zweiten Summanden gilt. Nun ist: 
Fri 63 | 
2 EIER nam n—eite—asnantm) | a 


t 


x —a\ (8 —o) pj2 
= [eenm nern unser nl | 
—0 
ut a 2% p .. 1: sin(#—0)pj4 
Fe SR a UNO Er (x —a) p]4 


Dieser Ausdruck hat ebenfalls bei &—= «a sein größtes Maximum, und zwar mit dem Wert = amp. 


Somit ist der Hauptvektor für alle x vernachlässigbar klein, die nicht in der Umgebung von 
x = a liegen. Nun besteht der Hauptvektor y(z,q) («+a) aus zwei Vektoren, die sich in 
Abhängigkeit von g beide mit der Winkelgeschwindigkeit —«a um den Nullpunkt drehen, wenn 
man sie in einer komplexen Ebene vom Nullpunkt aus aufzeichnet. Der Absolutbetrag des 


ersten Vektors ist (aq + b) Ta ;er ist also mit g veränderlich, während der Absolut- 
\ . (ea) pi2 7 1. sm(®&—o)pf4., | 
betrag des zweiten von qg unabhängig, nicht größer als Dr a ar ist und für —a 
den Wert = ap annimmt. 
Der resultierende Vektor dieser beiden Vektoren, der ja gleich dem Hauptvektor ist, wird 


also um so genauer für einen Wert x in der Umgebung von &= «a sich mit derselben Winkel- 


geschwindigkeit —& drehen, je weniger der erste Vektor sich im Verhältnis zum zweiten mit 
seiner Amplitude ändert. Nun variiert für x—=« die Amplitude des ersten Teilvektors von 
y(x,g) bei einer Verschiebung um den Betrag Ag um aAg. Diese Größe muß also klein sein im 


‚Verhältnis zu D ap, d.h. bei gegebenem a muß p groß sein, was wiederum nur möglich ist, 


wenn b genügend groß ist. Daraus folgt, daß man durch die Periodogrammanalyse einer Welle 
mit konstanter Frequenz und linearer Amplitudenänderung unter den angegebenen, ziemlich 
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Bild 2. Phasendiagramm von Kurve in Bild 1. Bild 3. Amplitudendiagramm von Kurve in Bild 1. 


allgemeinen Bedingungen die richtige Frequenz & aus der Neigung des Phasendiagramms der 
Versuchsfrequenzen x, in der Nähe von a bestimmen und damit nach Gl. (2) aus dem Ampli- 
tudendiagramm auch den Amplitudenverlauf berechnen kann. 


8*+ 


ine Gerade mit 


diagramm der Periodogrammanalyse der Kurve in Bild 1. 


Analysenintervall:p = 61 mm 5,3 mal enthalten ist, muß als Neigung im Pha 
erscheinen. Man mißt für n,,n;,n, der Reihe nach die Werte: 5,39; 5,35; 5,36. 


Abbildungen berücksichtigt. ar ZEN, Do 


Ein Beispiel möge die abgeleiteten Verhältnisse vera 


raden für die Versuchsfrequenzen: n—= 4,5 und 6. Da in Bild 1 die \ 


eine Übereinstimmung innerhalb 1%. Die Neigung im Amplitudendiagramm Bild st: 


mit der Neigung der Amplitudenfunktion in Bild 1, wenn man den Verzerrungsfaktor d 
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82. Die Periodogrammanalyse einer Kurve konstanter Frequenz “ 


De 


und exponentiell veränderlicher Amplitude. 3 


"Der Wellenzug habe die analytische Form: y = e-“tcos (at —P). Die verallgemeinerten “ = 


Fourierkoeffizienten in komplexer Schreibform sind dann: En \ 5 


.a+» a+» | Be: 


a(x,g)-+ib(z, =. | ertt te-nı- det nat+l | erertitetnites— md: u 
10 | p 


» ö ‘q : 2 
VIER 
. 2l-e+i@—a)] 
re op +ag—P] en] = y(x,9) +n (2,q) 


Im Hauptvektor y(x,g) steht in der geschweiften Klammer die Differenz zweier Vektoren, 
deren Absolutbeträge von q unabhängig sind und deren Argumente linear von g abhängen 
mit dem Proportionalitätsfaktor —o. Daraus folgt, daß das Phasendiagramm des Hauptvektors 
für jede Versuchsfrequenz zeine Gerade mit der Neigung —.a darstellt, d.h. die Frequenz kann 


. aus dem Phasendiagramm bestimmt werden für alle &, bei denen der Störvektor n (x, g) gegenüber 


dem Hauptvektor, absolut genommen, vernachlässigbar kleinist. Um dies zu entscheiden, bildet 
man den Quotienten d der Absolutbeträge beider Vektoren an der Stelle «=. Man erhält: 


Iyie q) a + er?eP —2 e-eP cos («—a) p} a (4 a), I 
Sr “Belo (2a) 

} ee 
:|n (®, q) ve au ES RA—A Erg cos (at) pt. re 


Eine Abschätzung lehrt, daß für d folgende Ungleichung gilt: 
i mer ? 
a ul Ne RR RE SE (5). 
a te | 
Setzt man hier « = : rund = — &, so erhält man für die rechte Seite: 
i RR ee 
ae Ye+4 "1 —_e-2ars 
Dieser Ausdruck strebt mit &—0O gegen en und mit e—>oo gegen den Wert 1. Die Aus- 
rechnung lehrt, daß für r>3 und e<0,2 der Wert d’<0,1 ist. Daräus folgt: Wenn das 


Analysenintervall » mindestens 3 Wellenlängen der Ausgangswelle umfaßt und e = En 0,2 ist, 
ES 


so ist der Störvektor praktisch vernachlässigbar klein. Die Einschränkung für & bedeutet, daß 
beim Fortschreiten auf der i-Achse um eine Wellenlänge die Amplitude der Welle nicht um 
mehr als 28% ihres Wertes fallen darf. Es bleibt noch zu zeigen, daß der Hauptvektor für 
solche Werte von x, die nicht dem Wert x benachbart sind, vernachlässigbar klein ist, gegen- 
über seiner Größe für die Werte von x in der Umgebung &. Durch Differentiation von Gl. (4a) 
zeigt man, daß der Absolutbetrag des Hauptvektors an der Stelle x —= & ein Maximum besitzt 
und alle übrigen Maxima symmetrisch zu dieser Stelle liegen. Die Tatsache, daß an der Stelle 
»=.a das Hauptmaximum liegt, ergibt sich leicht aus folgender Ungleichung: 
em2ea(] —e-e?)% , er?ea[| Her2er —2Derercos(e —a)p], e2ee(1 022)2 
Plot] PIET (@—a)] Sale 
Die beiden äußeren Formeln dieser Ungleichung stellen Kurven dar, die einander parallel im 


A er204. ee» 
Abstande (oberhalb der x-Achse) verlaufen. Jede dieser beiden Kurven be- 


[oe + (x —o)2] 


’ 
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arten B] era BER Ei 


Ei hau 


- wählt werden kann. 


ana) = ern lie} + nem. Ffnila + jeher N 
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Ra sitzt bekanntlich als Funktion vonxbeiie=xaein einziges Maximum und ist nur in der Nach . 
barschaft x = & wesentlich von Null verschieden. Der Abstand beider Parallelkurven wird 


für große Werte 0 p unendlich klein. Dies ist grundsätzlich erfüllbar, da p genügend groß ge- 


Zusammenfassend kann man feststellen, daß unter den angenommenen Bedingungen im 


.  Phasendiagramm für alle Versuchsfrequenzen x in der Nähe von x—= a gerade Linien mit der 
Neigung = —« entstehen, d.h. die Frequenz « ist aus dem Phasendiagramm bestimmbar. 


Das Amplitudendiagramm verläuft nach demselben Exponentialgesetz wie die Grundwelle. 
Nach Kenntnis von a sind also die Anfangsamplitude und der Exponentialabfall berechenbar. 


83. Theorie und Praxis der Periodogrammanalyse von Kurven 
konstanter Frequenz und beliebig veränderlicher Amplitude. 

Es seiz2= f(t)>0 der analytische Ausdruck der integrierbaren Amplitudenfunktion einer 
Welle mit der Kreisfrequenz & und der Anfangsphase ß; d.h. es sei: 


y= fit) cos (at -ß) BR LAN dt N ; ER: (7) 


a+p WR 
Rand + baden | sWec-mntten-nt snectmras . 
M q | | 


. Der erste Summand ist der Hauptvektor y (x, g) und der zweite der Störvektor nie, g). Aus- 


schlaggebend für die Lösung des vorliegenden Problems ist es zu wissen, unter welchen Vor- 
aussetzungen für f(f) die Gl. (8) folgendes unbestimmte Integral besitzt: Kar 


wiederholte partielle Integration zeigt: 


a+» i a+P 
1 REN A Le | RER USER NL 
z/10. ( rai=le ( Re ae p' 


- 


die Gleichung der Welle. Ihr verallgemeinerter Fourierkoeffizient lautet in komplexer 
Schreibweise: G Ne 


Dies ist offenbar erfüllt, wern f(t)>0 entweder monoton zu- oder abnehmend ist. Es gilt 
. aber auch für den Fall, daß f(t) eine ganze rationale Funktion n-ten Grades ist, da man durch 


Den Wert für 7(&,g) erhält man aus dieser Gleichung, indem man in ihr « durch —o ersetzt. 


Da sich praktisch jede aperiodische Amplitudenfunktion f (£) durch eine ganze rationale Funk- 
tion beliebig genau annähern läßt, so ist es tragbar, sich auf das Vorhandensein einer der drei 
genannten Fälle für das Aussehen von f(t) zu beschränken. Die Auswertung von GI. (9) liefert 
dann für den Hauptvektor 


ind = Flat ni a} et nF lgite—aljeitsetn. NR ON 


Für den Fall e=« muß der Hauptvektor existieren und sein Absolutbetrag ein absolutes 
Maximum besitzen, da folgende Ungleichung gilt: 


4+» 49 
frWee-nras| s/Fa)dt. 
g q 


.Aus Gl.(10) erhält man den Störvektor n (x, q), indem man in ihr den Wert «x durch —« und 
— ß durch + ß ersetzt. Daraus folgt sofort, daß der Absolutbetrag des Störvektors an der 
Stelle x = a kleiner sein muß als der des Hauptvektors. Dies bedeutet für das Phasendiagramm 


folgendes: Für = « nimmt Gl. (9) die Form an: fr 


md ra4+ mM —r@leretn + [Fat m 20) Fa 2a |ewe-n. 


Da im Argument des Hauptvektors die Variable g linear mit der Kreisfrequenz —a als Faktor 
vorkommt, so muß im Phasendiagramm für x = a die Kreisfrequenz — a als Neigung erscheinen. 
Denn, wie in der Arbeit Blumeund Kolb l.c. nachgewiesen ist, kommt im Phasendiagramm 
bei Überlagerung zweier Wellen auf derselben Harmonischen grundsätzlich diejenige Frequenz 
zur Anzeige, deren Vektor die größere Amplitude besitzt, und das ist hier, wie bewiesen, der 
Hauptvektor y (a, g). 

Um aber die Frequenz —«a aus dem Phasendiagramm allein erkennen zu können, ist es 
offenbar erforderlich, daß mindestens für alle Versuchsfrequenzen x, die in genügender Nähe von 


° 


x» = aliegen, ebenfalls die Neigung — a angezeigt wird. Dah 
gestellt werden, unter denen dies erfüllt ist. Zu dem Zwecke ‚der kt m 
vektoren, die nach Gl. (10) den Hauptvektor Y(mq) bilden, in einer 
Nullpunkt aus aufgetragen. A 


re 2 ‘von y(x,g) mit der Winkelgeschwindigkeit —o um.den Nullpunkt, da il 
ander von g unabhängig ist und q linear mit dem Faktor —« in ihrem . 


> Infolgedessen wird sich der resultierende Vektor y (x, q) allgemein nicht mit der Winkelges 
digkeit —« drehen. x 


Wellenlänge er ER. D= 2% und r groß, so wird für solche VersnensBeRBERBER az 
N " I « & « Ex 7. Id EN HE 
nen (m —=1,2..); die in der Nähe von « = oliegen, der Ausdruck (x —.4) = ER% if) n 


klein, und man kann setzen: 


Aus dieser Form ergibt sich aber folgendes: 


dabei konstant bleibt. 
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Beieiner Änderung von q um den Betrag A q > 0 drehen sich dann 


Im allgemeinen aber ändert sich dabei das Verhältnis der Absolutbeträge bei 


a TERRA 
Ist jedoch das Analysenintervall p groß gegenüber der zur Kreisfrequenz & gehö 
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Für derartige Versuchsfrequenzen kann also die Gl. (10) folgendermaßen geschrieben werden = 


innen rate) >.» (108) 


Bei einer Änderung von gum den Betrag A q wird sich y (x, q) mit der Winkelgeschwindigkeit —« 
um den Nullpunkt drehen, wenn der Ausdruck: Ya 


Ki Rat rg), ee. (10 b) 


Dieselbe Drehung wird erfolgen, wenn F(g-+ p) Sr 9 (4, p) [F* (p) —k], wo k 
eine Konstante und @ (g, p) eine Funktion von g mit dem Parameter p bedeutet. Da die Diffe- 
renz auf der linken Seite der Flächeninhalt zwischen f(t) und der t-Achse und den Ordinaten 
in qgund q-+ pist, so gilt der letztere Fall dann und nur dann, wenn f(f) eine Expontential- 
funktion ist. Er ist im $2 behandelt worden. WERT: 

Die Größe f(t„) nach Gl. (10b) stellt den arithmetischen Mittelwert der Amplitude der 
Amplitudenfunktion f(f£) über dem Intervall p im Bereich qg und g—+ p dar. Wenn f(t„) bei 
einer Verschiebung um Ag >0 also unverändert bleibt, so dreht sich % (2, g) mit der Winkel- 
geschwindigkeit —«. Für die Praxis ist es offenbar hinreichend, wenn f (t„) bei einer solchen 
Verschiebung sich relativ wenig ändert. Ein solcher Fall ist im ersten Paragraphen ausführlich 
besprochen worden, wo f(t) eine linear abfallende Funktion war. Der Fall einer relativ ge- 
ringen Anderung von f(t„) bei einer Verschiebung Ag kann aber nicht nur vorkommen, wenn 
/(t) eine monotone Funktion ist. Es kann z. B. sein, daß f (£) viele Schwankungsstellen im Bereich 
q bis q + p besitzt und dennoch f (tm) sich bei einer Verschiebung A gq< p relativ wenig ändert. 

Gerade dieser Fall kommt bei Registrierkurven aus Erschütterungsvorgängen häufig vor, 
wenn man das Analysenintervall genügend groß nimmt. (Siehe das Schlußbeispiel in der Arbeit 
Blume-Kolbl.c.) 

Damit ist es mathematisch klar geworden, warum man im Phasendiagramm auch dann 
eine verhältnismäßig genaue Frequenzanzeige erhält, wenn die Amplitude im einzelnen ziemlich ° 
schwankenden Verlauf hat. 

Im Amplitudendiagramm, wo f(t„) als Funktion von g angezeigt wird, kommen alle die- 
jenigen Eigenschaften der Amplitudenfunktion f(t) nicht zum Vorschein, die bei der Mittel- 
bildung,verlorengehen, d.h. aber die Gesetzmäßigkeit im großen zeigt sich auch in ihm. Ist 
z. B. f (t) eine ganze rationale Funktion n-ten Grades, so ist auch f (tm) eine solche; oder ist f (t) 
eine Exponentialfunktion, gilt dies auch für f(t„) usw. Man kann daher aus dem Verlauf des 
Amplitudendiagramms auf die Gesetzmäßigkeit von f(t) im großen zurückschließen. 

Ganz besonders vorteilhaft ist die Periodogrammanalyse, wenn es sich um das Problem 
der Abschirmung von Erschütterungen handelt, da hier eine möglichst genaue Kenntnis der 
Frequenzen erwünscht ist, während es genügt, vor den zugehörigen Amplituden ihre Mittel- 
werte zu kennen, um abschätzen zu können, ob eine Frequenz in ihrer Wirkung gegenüber 
einer anderen unberücksichtigt bleiben kann oder nicht. 

Der Fall f(t) = konstans, für den allein die Gesetzmäßigkeit des Verhaltens im Amplituden- 
und Phasendiagramm abgeleitet wird, ist einfachster Sonderfall der angegebenen Theorie. 
Eingegangen am 7.2.44, 


Eine wahrscheinlichkitsheoreische Behandlung dı der 
EL | |  Jokereigenschaft. Ey Ne 
BERN 2 a y Be Von Hans Münzmer i in Göttingen‘). 2 B SE 


BE, Die. AN daß in einem ‚Spiel einer Karte die Gültigkeit einer beliebigen. Karla Pe ( Join Prag 
E. 4 R m ‚Romme), wird einem Wahrscheinlichkeitsmodell zugrunde gelegt und die Wirksamkeit dieser Bigen- 
‚schaft an einem Beispiel gezeigt. Die Ergebnisse werden Buch das Veen u Kundenwerbung durch de 
. Packungen beigelegte Gutscheine angewendet. 


Een: Q Sr The property that the value of one card of the ar is to TA Hat: of any tarlio one 8 Hiking (e. g: the 
re, ““joker‘‘ for Romme) is taken as basis of & model of ‚probabilities and the, SAY S this property ı Lee 
ii demonstrated by virtue of an example. R 


The results are applied to the practice of puffing ome’s goods by extra tickets enclosed in 'packages. p 


a ” 
ey... f 17 propriete, pas laquelle une ceriaine carte est succeptible de ‚prendre la valeur de n’importe quelle 
re  auire carte du jeu (p. exp. »joker« dans poker) a &t& prise pour base d’un modöle de probabilite et son effi- 
cacite illustree Ber ‚un ewemple. Des neoullaia sout SER a une methode de reclame Bar les bons joints 
aux paquelts. i } 
B ocHoBy Mozenn Bepoarnocru NONOKEHO TO CBOHCTRO kaprounoN urpit, "KOTN& onpe- 
‚  AeleHHaA Kapra (HANp. ‚‚mKoKKep‘‘ B TOKepe) MOkKeT IIPHHUMATL 3HayeHNe moboänsuroh 
‚. KAPTbl KONONE, UpHYeM EhheRTHBHOCTL HTOTO CBOÄCTBA TOACHeHA HA mpumepe. Pesyasraru 
HpPHMeHeHBI K MeTOAy PeKHaMEH UYTeM BEIAAEIBAHUA BEIUTPBIIIHEIX TAJIOHOB B JHaKOBKY. 


Jakob Bernoulli erwähnt im 18. Abschnitt des 3. Teiles seiner: ‚‚Ars conjeetandi® 
(1713) verschiedene Kartenspiele, bei denen einzelne bestimmte Karten nach Belieben an Stelle 
jeder anderen Karte des Spieles verwendet werden dürfen. Es werden das deutsche Treschak (auch 

. unter Bretling oder Krimpelspiel bekannt), das französische Brelan und das spanische L‘'hombre- aa 
Spiel genannt. Gegenwärtig ist von derartigen Kartenspielen am meisten das Romme verbreitet, An 
bei dem den Jokern die Gültigkeit einer jeden beliebigen Karte zukommt. Wir werden diese 
Eigenschaft, die wir als Jokereigenschaft bezeichnen, einem einfachen Wahrscheinlichkeitsmodel 
(Urnenschema) zugrunde legen und dann ein Problem behandeln, das besonders deutlich die Wirk- 
samkeit dieser Eigenschaft zeigen wird. In einem der Beispiele wird dann das gewonnene Ergebnis 
auf ein Verfahren der, Kundenwerbung, nämlich auf die Werbung durch Gutscheine, die un 
Warenpackungen beigefügt werden, angewendet. 

E; Eine Urne enthalte a mit einer Null und je m mit den Nummern 152, ‚n beschriftete Ku. 

Be ' geln. Ihre Gesamtzahl sei mit N = nm + a bezeichnet. Die Ziehungen, die a aus der gut durchge- 

_ mischten Urne vorzunehmen sind, sollen ohne Zurücklegen der bereits gezogenen KuRen erfolgen, 
so daß Abhängigkeit zwischen den einzelnen Ziehungen besteht. 

Wir fragen nach der mittleren Anzahl H (n, m, a) = H der Ziehungen, die man machen muß, 
bis man eine vollständige Serie der Nummern 1,2, ..., n erhalten hat. Dabei soll gestattet sein, die 
gezogenen Kugeln mit der Beschriftung Null, die wir im folgenden kurz Nullkugeln nennen wollen, 
an Stelle beliebiger noch fehlender Nummern zu verwenden (Jokereigenschaft). Es liegt also z.B. Beh: 
bei n = 50 schon mit 46 verschiedenen Nummern und 4 Nullkugeln eine vollständige Serie aller 3 
Nummern vor. Außer dem Mittelwert H wird noch der zugehörige mittlere Fehler o ermittelt. Da 1. 
es sich hier um die Zufallsveränderlichkeit der Ziehungsanzahl s bei festgelegtem Ziehungsergebnis Be 
(eine vollständige Serie) handelt, gehört das hier betrachtete Problem dem Pascalschen Pro- | 
blemkreis an. 

Die folgende Berechnung besteht im wesentlichen in ‚der Bestimmung der Währseheinliehkeit 

- w (s), daß in s Ziehungen aus der Urne keine vollständige Serie der Nummern 1 bis n erscheint. Für 
s <nistw(s)=1. Die Differenz w (s —1) —w (s) ist dann die Wahrscheinlichkeit P (s) dafür, 
daß gerade mit der s-ten Ziehung die Serie vollständig wird. Der gesuchte AU aUeE ergibt sich 
bei Beachtung von w (N) =0 auf Kate der en Definition als 


Hnm,a)= SP (= Zstwe— —w(]= Luis FSB A RR ZA NP 


vr 


N a un deh  d Ze 


Für das mittlere Fehlerquadrat o? erhalten wir 
’N N N 

= 5 (s—H)2 Pi) Fs2[w(s—1) —w(s)] — = 28 (s+1)w() —H—H! . . (2). 
s=0 sl s=0 


Die Wahrscheinlichkeit w (s) führen wir zurück auf die leichter anzugebende Wahrschein- 
lichkeit % (s, k, r), daß in s Ziehungen gerade k Nullkugeln erscheinen und dabei r verschiedene, 


bestimmte Nummern fehlen. Von —— at .) möglichen Fällen sind hierfür Aus 


re) Fälle günstig, so daß wir 


1) Aus dem Institut für mathematische Statistik und Wirtschaftsmathematik der Universität Göttingen. 


aa Mit v ER Mb er wir nun nie Walrsche li | aß in 
kugeln gezogen werden Ins mindestens k+ 1 verschiedene Nummern fe 
Ban etabh, sich n ER 


ER RE N RE N Re) r—1 PERS 2 wi u IR 
Re IR “ x eis name = ( e) ehe) en N: AR 2 


"Das erste Glied er rechten Seite von (4) stellt die Summe der Wahrscheinlichkeiten dar, 


das Fehlen der n: )v verschiedenen &-+ 1-Tupel bestehen. Es lautet daher (ei “ 1)* (8, hu | 
x Dabei ist für ‚jedes, der. 2: en Me verschiedenen k+ 22 Tupel die. Wahrscheinlichkeit & 


u(s,k,k+ 9) k A 2-mal mitgezählt worden, muß also k + 1-mal wieder abgezogen werden, da sie 
 inv(s,k) nur einmal berücksichtigt werden darf. Damit ist das zweite Glied der Formel (4) erklärt. 


Die Wahrscheinlichkeit u (5 k, k.+ 3) ist nun für jedes k+ 3-Tupel i im ersten Glied (3 sr >)- mal 


“ mitgezählt undi im zweiten Glied 7 > ) u * >) -mal abgezogen worden. Dassie in v Ve k) für jedes - ‘ : 


der k u. 3 verschiedenen k zu 3-Tupel wiederum genau einmal eithalten sein muß, ist als anittes ; 


ai ELEH-CHN slemanın- (et )@t euere 


2 azaardeen. Durch dieses abwechselnde Zuzählen- und | Ahtichen entsteht Formel (4). en ist für a 
0<k<n—1 gültig und geht für k= 0 in den bekannten Additionssatz Seh der Wahr- N 
scheinlichkeitsrechnung für sich gegenseitig nicht ausschließende Ereignisse über. u 

Da v(s, k) als Wahrscheinlichkeit dafür aufgefaßt werden kann, daß in s Zehen kNul- 
Kugeln erscheinen und dabei keine vollständige Serie aller n Nummern zustande kommt, ergibt sich 
die totale Wahrscheinlichkeit für das NICHLENS ENGE un einer Bene See in s Ziehun- 
gen als 


E | | | (8) = ARE BSR RAN nn ie; 


Bei Beachtung von (9) — 1und | I —0fürd<0 und b> agilt diese Formel für 0<s = N, 


so daß (1) ohne Einschränkung angewendet werden kann. Setzen wir (3) in (4), (4) in (5) und (5) 
in (1) ein, so erhalten wir den gesuchten Mittelwert Sn ee! 


Eat 


NEE 6). 
s=V k=V r=k-+1 ER . 
1; g 
a Die Summation über s läßt sich hierin, wie folgt, durchführen: | 
Br: j e N ar | 
f ke nm —rm ;) re nm+ta+l ! 1 
Et. g| 8s—k ) Ze Er BER Ba oo nm+a+1 7) 
Be —, (nm+a rm-+a\ [nm+a\ /(rm-+a\ (nm ; 
) ee IE u 


wobei die erste Gleichsetzung nach Darstellung der Binomialkoeffizienten durch Fakultäten sofort 
einzusehen ist und die Summierung selbst nach einer bekannten Formel der Kombinatorik (vgl. [3] 
Seite 251, Formel (26)) )) erfolgt. Y 

a wir (7) in (6) und ersetzen wir die Summationsgrenzen k+ 1 <r<n und 
0<k<n —L-durch die neuen Grenzen 1<r<'n und 0 <k<r —1, was entsprechend dem 


Dirichletschen Satz über die Vertauschung der Integrationsgrenzen bei Doppelintegralen er- 
laubt ist, so erhalten wir 


| t Ka: r—1 | 
a Je 2 8) 


| m +a EN 


j EM NR N BB h Ne) er Fe en i | ne | 
un En ba Pan | B "mai = rs 5) nun auch die Summation über all 

SEN En NEN u Re 3 k H 2 L ar ’ e + Br j AN ke 
Dan vorgenommen werden (vgl. BL: ‚Seite 252, ‚Formel e)). Hg / - 
EI Pr Es reiht sich schließlich * | Ps; TE 

ren N ya a 

en CH n, m, 1 e> mt mr HEN] ee 

Et ‚En L Ma) — ae BL -1) AR | mtart a Re 

” Wir: 2 | e er | a ii, ur ' ale (9). 

ae mtr) Fa ER EINEN 


nr ee entsteht aus 8) (4 4,6 (5 )und 2) das mittlere Fehlerquadrat, das unter Ver- K En 
Aa wendung der Hilfsgrößen ! an, 


N Mer Han N ER Ri | ner ; 5 Ss & a 2 
Br, ? ng u % NEN, #% we 


"sich folgendermaßen darstellen läßt: DS VER: 


F- Ber B =20mtat mat Ic ten —H(n,m,a a) —H?(n,m,a) (11). 
= i el 


"Von Interesse ist der Grenzübergang Mm—> ©, a— © bei festgehaltenem Bee Nach kurzer 
Rechnung erhalten wir aus (9) _ 4 Ka 
Be af | 
4 a lahm, m, mb) =H (n, b) I 1)r . a ER (12) , 
und ebenso aus (11) ri BE 
ie —2 
; Im Sylt —_. EZ: > re b)—H( mb) ZHen,. Sn 3). 
ER 
+ 


- Man kann diese Grenzformeln auch direkt herleiten, wenn man von DR Urne ausgeht, bei 


\ M>. 


der die für alle Ziehungen unverändert geltenden Wahrscheinlichkeiten p = + ? für das Ziehen 
einer Nullkugel und qg = - für das Ziehen jeder der Nummern 1, 2, ...,n bestehen. Die 


= Ziehungensindhierim Gegensatzzumzuerstbehandelten Urnenschemavoneinanderunabhängig,wass 

praktisch durch Zurücklegen der gezogenen Kugel nach jedem Zug erreicht werden kann. Die direkte Bo 
Herleitung verläuft genau so wie imallgemeinen Fall, lediglich die hypergeometrische Wahrschein- Br, 
lichkeitsverteilung (3) ist durch die Binomialverteilung zu ersetzen. 

Durch Umformung erhalten wir aus (9) und (11) beim +0 für H (n, m, a) und o? andere Aus- | 
drücke, die einerseits den Anschluß an die bereits bekannten Formeln des Spezialfalls a = 0, m— © 
gestatten und andererseits für die numerische Berechnung bei größerem r und kleinerem a von Vor- . rL 
teilsind, da sich die Summationen nunmehr über « + 2 bzw. a + 3 Glieder statt über n Glieder er- 
strecken. Wir geben hier gleich die Endformeln an. 


Nr “ % 
IF En iu 


a+1 


> Be aD Cole; BEER: (14) 
a - m—1 m—1 m—1 
al er ) & (" nm) IR T5) 
een 
un 
a Sc: 1)» Pr a rer ) E-—H2.(16), 
2 en 


Bar 
a 


wir die einfache Formel" 
Er (n,m, DE 


ji 


% 


übergeht CH) [2], LA. Das dazugehörige mittlere Fohlerquadrat Tantet dann für 
| en Er 


ee) 


N vehohni en Urmformungen aus a 6) u Für den Grenzwert ergibt RL wie bereits be 
ö Rau, B t ne Ne ren BE RN, 


| r Sr E N 
RM, E 


RER im ot— nt ) ar — Re ee NEL 
4 1: Beis pie 1. Wieviele Rarlen muß man im Mittel einem Haie Iurcngeniienten Quartett- ee 
spiel von 12 Quartetten entnehmen, bis man von jedem Quartett mindestens eine Kartehat? Da- 
 beisollen noch 0 Joker, bzw. 2 Joker, bzw. 4 Joker zunatzlich im ‚Spiele sein. Na den Formeln Te 
ebis 1) oder (14) bis (16) erhalten wir | } IRRE a 


DR "2 (1,40) = 35,5 mit, Rest | 
FRA 1.282,49) >205 mito— LUA, 
ei E24 - 195 mit = HEN 


2. Beispiel. In der Kundenwerbung hat ein Verfahren immer mehr Nase ge- Ve 
‚ funden, das das Spiel des Zufalls in reizvoller Weise benützt. Den Warenpackungen zahlreicher Ge- 

brauchswaren werden Bilder oder Gutscheine beigelegt, die die Käufer sammeln sollen. Beidem 

allgemein beliebten Sammeln der vor dem Kriege den meisten deutschen Zigarettenpackungen bei- 

gegebenen mit Nummern versehenen Bilderschecks handelt es sich darum, vollständige Serien der har 

Nummern 1 bis 50 zu erlangen, Der Anreiz zum Sammeln wurde durch sogenannte Freischecks ' = % 

erhöht, die in einem gewissen Prozentsatz vorhanden waren und entsprechend den Nullkugeln 

unseres Urnenschemas jede der 50 Schecknummern vertreten konnten. 

Nimmt man an, daß 1. P%, der Gesamtzahl aller Schecks Freischeck waren, 2. die Nummern- 
schecks mit den Nummern 1 bis 50 gleichhäufig vorkamen, 3. die Gesamtzahl aller Schecks sehr 


0... groß war und damit praktisch konstante Wahrscheinlichkeiten für die einzelnen Schecks vorlagen, 
0,4. die Verteilung der Schecks auf die Packungen ganz nach dem Zufall erfolgt war, 5. schließlich 
See kein Tausch der übriggebliebenen Schecks zwischen den einzelnen Sammlern stattfand, so gibt die 
je. Formel (12) die mittlere Anzahl der Schecks an, die man bis zum Erhalt einer vollständigen Serie 
I RN a U als 
BR nn ee Ki 
= = | - | Be erwerben mußte. Dabei ist b= T00—P A Beh, 
SEE RDETETET PERRTEUREITRT TEFAL" NURKHT ARE Aus Die nebenstehende Tabelle enthält auf ganze Zahlen 
 ; 0 0,0000 225 463 abgerundet für verschiedene Prozentsätze diese Mittelwerte 
@ 1 0,5051 184 + 45 mit den dazugehörigen mittleren Fehlern. Sie zeigt, in ° 
8; 2 1,0204 160 + 33 ‚welchem starken Maße der Mittelwert 225 durch Beigabe 
N ne Ser <ü es von Freischecks vermindert wird. 
20 12.5000 778 Wenn gleichzeitig mehrere Serien gesammelt werden, 
50 50,0000 DIEB kommt ein beachtlicher Ausgleich zustande, wie für P=0 
100 @ 50+ 0 in einer früheren Arbeit [2] gezeigt worden ist. 
Literatur. 


[1] H. Münzer: Über eine einfache Lösung einer Wahrscheinlichkeitsaufgabe in der Kundenwerbung. Z. 
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Eingegangen: Mai 1947, 


| B koordinaten?) 


N 2 KTEINE MITTEILUNGEN 


Die ande dünner Kreisringplatten, die am 
Außen- und Innenrand gleichmäßigen Drück 


‚erfahren 1), 


Wird eine ones der Dicke h am eh 


Y and Innenrand durch einen überall mi Druck p 


beansprucht, so hat man als Platteng eichung in Polar- 


4% 


| kur 0, 


wo w die Ausbeulung aus der Plattenebene mißt und 
el (F 


Kr 


"Ram 


Ei z. „Plattenfaktor ‘) ist. - Die 


| Randbedingungen lauten dabei 


a) beieingespanntem Rande 
| iw_g 


ee 
2) dr 
b) beidrehbargelagertem Rande 
dw “ 
ZENARRE el =,00.. (2b). 
w-0, Ze ne (2b). 


“ Führt man y = 2 als neue Veränderliche ein, inte- 
r 


 griert einmal und dividiert dann noch durch r, SO er- 


‚hält man aus Gl. ia 


d | C 
Anus EU 
a 
Man sieht leicht, daß 
a ee 4 
M: Hr (4) 


3 ein pertikuläres Integral dieser Differentialgleichung 
ist. Nach Multiplikation mit r? geht Gl. (3) durch 


‘die Substitution z—= YA: r, c,— c,|VA über in 
„0° dy 
ee „te-be=az.- 


eine Besselsche Differentialgleichung. Da iı in Gl.(4) 
schon ein partikuläres Integral gegeben ist, genügt es, 
die homogene Gleichung zu integrieren. Man erhält” 


va=ah@d+GN,@ .... (6) 
Das allgemeine Integral von Gl. (3) lautet daher 


5) 


—— c 
vn=GHhNMn)+C,N, (Var) + 3% . 
und daraus ergibt sich durch Integration 


wer)=k,) VAr)+k,N(YAr)+k,inr+%, (8), 


Für die Randbedingungen nach Gl. (2a, b) am äußeren 
(r =r,) und am inneren (r—=r;) Rand werden noch 
die folgenden Ausdrücke benötigt: ' 


d 1 kz 

1 5 —|r, VAJ, (VAr) + kz YAN, Vin] (9a) 
drd 

u a nu|- | 

— la KR A—T EnN 


+ kalarN. (Var) — VA —u) N, VAn)] 
A u EEE, 


(7), 


Im; Falle eines freien Randes ist hz =. 


1) Mitteilung des Lehrstuhls für angewand te Mathematik (Prof. 
Dr. Willers)an der Techn. Hochschule Dresden. 

») Fr. A. Willers: Die Stabilität von Kreisringplatten. 2. 
angew. Math. Mech. 23 (1943), 8. 253, Gl. (4). 


berechnet: De} 


jlefert mit rifra= o und YAra=® 


meld (eo) No) — 


. (2a); 


»IH. Innenrand eingespannt, 


IV. Innen- und Außenrand drehbar. 
‚gelagert; j 
V.Innenrand eingespannt, Außen- 

rand frei; j 

VI. Innenrand drehbar 


"Aus Gl. (0) und (08; b). wurde a für folgende Fälle 


I. Innen- RE eingespannt. 


Aus den Randbedingungen erhält man ein Gleichungs- 
system für die k,...%y. Nullsetzen der an. 


mm un 
+ ooln old, (e o)N, (0) — I, (0) Nı(oo)] | 
FROM HL 


Für 0>0 erhält man durch Grenzübergang unter 
Her, der Hankelschen Formel?) 


4 2 +N (a ne (@) + 7,(0) (v+1n 3- N i 


le 
wo y= U ‚57721506. Br 
ist. | 
Ferner ergibt sich hier durch längere Rechnung 
| 2 BarR» 2 
En NEE RR TI ER ee LO'b) 0 u 
ae Fand, ur 3 
Zu Bacon‘ etwas Kunlnre Gleichungen, als 
es Gl. (10) ist, führen die weiterhin untersuchten Fälle 


I. Innenrand drehbar 
Außenrandeingespannt; 


Ei Euler sche "Konstante ; 


gelagert, 


Kos 
rand drehbar gelagert; 


gelagert, 

Außenrand frei. a i 
\ 

Setzt man im letzten Fall o®@= (l1—e)w und ent- 
wickelt die in der Gl. (10) entsprechenden Gleichung 
auftretenden Besselschen und Neumannschen. 
Funktionen in der Umgebung von nach Taylor, 
so erhält man unter Beachtung der Beziehungen 


n-—h A--h-h Ih 
N=—N, NE NENNEN, und Br 
: ir 
2 ed! 
J, (%) N, (&) —J, (X) N, (&) ee E%- 
T% “ 
in der Grenze für e>0 | 
a Eee an). Me. 


Es ist der einzige der oben behandelten Fälle, bei dem £ 
für o — lein endlicher Grenzwert existiert. ' 


Für die übrigen Fälle kann man aus den Formeln 
für die Beullasten eines unendlich langen Platten- 
streifens®) asymptotische Näherungsformeln für ® 
ableiten, die z. T. schon für kleine Werte von og von 
den wahren Werten nur um, wenige Prozent abweichen. 

Die numerischen Rechnungen wurden für den Fall 
u = 0,3 durchgeführt und ergaben die in der folgen- 


“ den Tabelle angegebenen, in nebenstehender Figur 


dargestellten Resultate. In ihr ist ® als Funktion 
von o aufgetragen. Eine zweite Skala gestattet, 
o®—Ar, direkt abzulesen. 


°s)G.N. Watson: A Treatise on the Theorie of Bessel Func- 
tions. Cambridge 1922, $. 62, Formel 2 spezialisiert auf n = 0 
bzw.n=1. — \gl. auch Math. Tables VI, Bessel Funetions, 
PartI: Functions of order zero and unity. Cambridge 1937, 8.174. 

4% Fr. A. Willers: Das Ausbeulen von Plattenstreifen, deren 
Dicke sich sprungweise ändert. Forschung auf dem Gebiete des 
Ingenieurwesens 10 (1939), 8. 231—232. 


6,648 | 6,648| 4,401 | 4,401| 2,049 | 


Re nr 
SU A| Ba | 


0, 6,815 | 6,454| 4,535 | 4,278) — 
0,1 TR ‚115 6,307 | 4,762 | 4,205 
2 ? 7,926  6,476| 5,360 | 4,341 1,801 
v3 ‚9,021 | 7,058|, 6,156 | 4,749 1,613 
0,4 10, 501| 8,001] 7,226| 5,410| 2 1,452 
0,5 112,586) 9,416| 8,726) 6,401 1,323 
0,6 .|15,721| 11,604 | 10,976 | 7,933 | . 1,220 
O4 | 20,953 15,305 | 14,725 | 10,523 1,135 | 
0.8 |s1.421 22,759 | 22,217 | 15,738 | 7,780 | 1,06 
0,9 | 62,833 | 45,199 | 44,688 | 31, ‚429 1 605 1,00 
Fe = & \ ) 


al 


2 


Art. 


a ı 
GER 94 06 Lr-% A 3 


Der Einfluß, den u auf ® hat, ist im allgemeinen 
sehr gering; z. B. ergab sich im Falle III bei o = 0,5 
die in der folgenden Tabelle gezeigte Abhängigkeit. 
Es muß jedoch berücksichtigt werden, daß u auch in 
den Plattenfaktor eingeht, und dieser Einfluß auf die 
Beullast P (= p-h) ist beträchtlich, wie die dritte Zeile 
der Tabelle zeigt. 


m 08.%0.1°°:0,2 ,7°0.32%0.4.220:5 
@ 8,606 8,647 8,687 8,726 8,765 8,803 


br Fe) 

— = _ 6,172 6,294 6,551 6,973 7,622 8,610 

Er 12(1—u%) | 
Sehr auffällig ist der Verlauf der Kurven II, IV und 

VI. Durch Verallgemeinerung eines Satzes von Cou- 

rant?°) ergibt sich nämlich folgender Satz: 
Sind die Eigenwerte Courantsche Maximum- 
Minimum-Werte eines quadratischen Variations- 
problems n-ter Ordnung und sind auf einem Teil 
des Randes unter den Randbedingungen diefolgenden 
enthalten: 


y=y-yr- 


so werden sämtliche Eigenwerte größer, wenn man 


an un) 2 0, 


5) R. Courant und D. Hilbert: Methoden der mathe- 
matischen Physik I. 2. Auflage, Berlin 1931, 8. 355. 


P den Fällen II, I 
Rand die Bedingı 


dem Gebiet, die Beullast. > ist ar 
aus plausibel, da mit kleiner werdender 


h | * N 
S ; h 03% 
1 B | 


Euer! 


mit kleiner werden 
Fall VI schon von Es — ST a A \ 


die Pastensteiigkei zunimmt. 
Dresden. 


EA/ Schuber 
Sun: 


N Eine Bemerkufe zur Neigungsempfindlichkeit 
.des Prandtlrohres bei Messung des stati schen ER 
Druckes. f 


Ein Prandtlrohr werde mit der nl 
keit vd. unter dem Winkel & angeblasen. Wir zer- 


legen d. in eine Komponente v, in Richtung dr 


Rohrachse (v; = v„ cos a) und in eine Kompo- 


Bons dv, senkrecht dazu (v, = v. sin 0; vgl. Big 2: | j 


 Druckmehring: 


D, : e e 
Bild1. Bild 2. 


Betrachten wir die Ehena des Druckmfinges a 
Bild 2), so ergibt sich als Betrag der Geschwindigkeit 


am Druckmeßring va —=2v, |sin@| und damit 


v= y» AERDe 
= Yo% costa + 4% sin? asin?-p oc) 
— Vo Veos®a + 4 sin? a-sin?@ 


Nach Bernoulli folgt daraus 


ee REN 
Q OD oo4) a a hen (2). 
MEZ 
(1) und (2) liefern 
PP —gin2g-(1—4sin?o) . . (B): 
2% | 


Integration über den gesamten Druckmeßring und 
Division durch 2% führt zu dem Wert, den das 
Prandtlrohr theoretisch anzeigen müßtet): 


+ 2n 
P—Po _ _ sin? a 


ee . (4). 


ER 
2 


Wir wollen nun zeigen, daß der gleiche Wert wie bei 


(4) auch in den folgenden Fällen theoretisch angezeigt 
werden müßte: 


I. Die Druckmeßstellen bilden einen kontinuier- 
lichen Halbkreis. Dabei ist die Lage der Druck- 
meßstellen zur Richtung von p, gleichgültig. 


Dieser Fall ergibt sich sofort aus 
S eg+n 
Zei SW —4sin?o)dp= —sin?2«: 
2 


!) Wir setzen hierbei voraus,daß sich der arithmetische Mittel- 
wert des Druckes einstellt. 


J L. D » 
ER 


E 
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II. Die Druckmeßstellen bestehen aus » (> 2) ein- 


zelnen Bohrungen, die gegeneinander um den ° 


gleichen Winkel “® versetzt sind. Auch hier 
kommt es nicht auf die Lage der Bohrungen 


hinsichtlich der Richtung von v, an. 


Um dies zu beweisen, betrachten wir auf Grund von 
(3) folgende Summe 


n—1 


Pro _ 1 ne 2R\|. sine 
ee Br 2 1 sin (2°) sin? a 
2.7 (8). 
_sin?o et 
4 aus +)| 
v=0 


Fürn =1 und.n = 2 ergibt (5) 


P —P® — sin? 

fe 

also einen von p abhängigen Wert, 
Zur Berechnung von 


a (1—4sin?p), 


n—1 
27 
> sin? (9 +) 
„‚=0 


müssen wir diese Summe in geeigneter Weise um- 
formen: 


n—1 


27 
Dr sin (++) 
‚=0 
n—1 


R 27 ß 2m\: 
= > sinpcosy + cospsinv 
‚=0 
"I in? 2 an 2 in? an 
= > sin? p cos vv, te psinvT, 
’=0 (6). 
i RTL In 
+ 2sinpcospsinv — co8y 


n—1 nl 
= sin? > Co + 008° p > sinay 2 
v0 0 
NR! 
2 2n 
+ sin29 > sinv — cosV Di 
v»=(0 
Wegen 
N 
sinvT, = — sin (n —») 5 
und 
FR 
c08y — = C08 (n — v) e 
ist die letzte Summe von (6) gleich Null. Zur Be- 
rechnung von 
eh N 
IR cos2y2 7 und Nam, 28 
N au! 
v=0 v=0 


gehen wir von der Gleichung‘ 
an —1=0 


aus, die die Wurzeln 


27 DIR 2n 
er — cosy-" +isinv ZT 
N 
n n 
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besitzt. vw =0,1,2, ..., #—1). Die Koeffizienten 
der Gl. (7) sind bis.auf das Vorzeichen die symmetri- 
schen Grundfunktionen ihrer Wurzeln, d.h. für n> 2 


n—1 


4v=0 (8). 


Nach den Newtonschen Formeln über den Zu- 
sammenhang zwischen symmetrischen Grundfunk- 
tionen und Potenzsummen folgt nun sofort 


n—1 n—1 2 
2n 
> von — ( > cos 
N n 
v=0 v=0 
n—1 
+2 5 EL eg; 
n n 
N) 
w<v 9) 
n—1 n—1 ( 5 
I" sine» 2 —( > sinv2) 
— N n 
»=0 De) 
n—1 
+2 2 sin we" r siny = — 0 
Av=0 
u <v 


Durch Subtraktion der zweiten Gleichung von (9) von 
der ersten erhalten wir unter Berücksichtigung von (8) 


n—1 9 3 
Sr ana) 
Su es sin en ur! (10). 
De) 
Damit ergibt sich 
n—1 2 
N" 0082» — > sind 277 
v=0 & „—— 
n—1 
en 2 
- = 
1 2 
und 
n—1 n—1 
N’ oosy2®= N siny2a—r 11 
Zr in? y 5 (11). 
v=0( v=0 
Dies liefert, in (6) eingesetzt, 
n—1 9 > 
7 
DL sin (o+,%)-3 . (12). 


»=60 


2 K Bi! 


anschaulicht diesen Sachverhalt: 


(5) und (12) führen schließlich m 0° 


PIALSN Fr 's BR ; l { 
O2 -Utm—2n]=—eina . .(13). 
% yo 2 \ | EEE 
N » r j F j 
chen a ars ihebrekiack Rt 
‚Vergleichen wir die theoretischen Werte von DEE 
”% 2, 


“ mit gemessenen?), so sehen. wir, daß die Meßpunkte 


von x—=0° bis « =20° sehr gut mit der Theorie 


übereinstimmen, bei höheren Werten von aber 


merklich abweichen. Die folgende Zahlentafel ver- 


Die Unterschiede bei größerer Schräganblasung können 
so erklärt werden; betrachten wir die Ebene des Druck- 
meßringes, so wird in ihr ein Kreis mit der Anström- 
geschwindigkeit v, umströmt. Nun ist bekannt, daß 
bei der Umströmung eines Zylinders die beobachtete 
Druckverteilung mit der theoretischen nur auf der 


- Vorderseite übereinstimmt?). Die ' Rückseite wird 


durch. Ablösungserscheinungen gestört. Wir können 
vermuten, daß dies auch der Grund für die Unter- 
schiede zwischen Theorie und, Messung in der vor- 
liegenden Zahlentafel ist. Der Einfluß dieser Ab- 
lösungserscheinungen wächst mit v,, d.h. mit dem 
Winkel &. Da nach I theoretisch bereits ein ‚‚halber. 
Druckmeßring‘‘ ausreicht, um den gleichen Wert des 
statischen Druckes wie ein gewöhnliches Prandtl- 
rohr mit ‚ganzem Druckmeßring‘‘ zu liefern, wäre 
folgendes Experiment aufschlußreich: man messe den 
Einfluß der Schräganblasung bei einem Prandtl- 
rohr, dessen statische Druckmeßstellen einen Halb- 
kreis bilden, und drehe das Prandtlrohr so um 
seine Achse, daß dieser Druckmeßhalbkreis einmal an 
der der Strömung zugewandten, das andere Mal an 
der der Strömung abgewandten Seite gelegen ist. 


Oberammergau. Hans-Joachim Kanold. 


2) Vgl. Wien-Harms, Leipzig 1932, IV, 1, 8.495, Fig. 11. Nei- 
gungsempfindlichkeit des Prandtirohres bei Messung des 
statischen Druckes mit scharfer Schlitzkante. 

®2) Vgl. Glauert-Holl: Die Grundlagen der Tragflügel- 
und Luftschrauben-Theorie. Berlin 1929, S.29, Abb. 14. 


Die günstigste Verteilungsbreite, ein neues . 


Streuungsmaß. 


Im folgenden soll ein Verfahren erläutert werden; 
welches bei der Auswertung statistischer Verteilun- 
gen verwendet werden kann. Insbesondere soll dieser 
Aufsatz zur Bekanntmachung eines neuartigen Streu- 
ungsmaßes dienen. Die Notwendigkeit der Ein- 
führung dieses neuen Maßes ergab sich aus praktischen 
Bedürfmissen. Seine Anwendung ist besonders für die 
schiefen Verteilungen gedacht. 

Neben einem speziellen Mittelwert (arithmetisches 
Mittel, häufigster Wert, Zentralwert), dessen Wahl 
sich je nach dem Zweck richtet, der mit der statisti- 
schen Erhebung verfolgt wird, läßt sich bei der schiefen 
Verteilung mit Vorteil der neue Streuungfaktor als 
Maß für die Anordnung der einzelnen Meßwerte inner- 
halb der Verteilung verwenden. Wie der wahrschein- 
liche Fehler beider normalen G au ss schen Verteilung, 
soll sich auch hier dieser Wert auf die Hälfte der Fläche 
unterhalb der Verteilungskurve beziehen, und zwar 


 Abszissenachse ermittelt w 


soll das Flä henst 
 szissendifferenz d 
. Fläche möglichst k 


also die Bern: 50% 


renz zu bestimmen. 


also. darauf hinaus, das Minimum d 
nv e 


e 


} ? 


BR 

Zu diesem Zwecke sei folgende analytische Btrac- 
tung eingeschaltet. Nebenstehende Skizze (Bild1) } 
stellt‘ die Integralkurve y=f(x) der Verteilungs- 


. funktion y’=f’(x) dar. 


x und z sind die Begrenzungsabszissen des Flächen- 
stücks, . ; 
q ist die sogenannte Verteilungsbreite, 

© = const. = 50%, ist der Ausdruck für die halbe 
Fläche. RE 

. Die Tangenten T, und T,in den P, und P, mit den 
Differentialquotienten y, und %’, veranschaulichen 
die Begrenzungsordinaten in der Verteilungsfunktion. 
Die Abszissendifferenz g ist definiert als g=z— a. 


Hier ist ; 
:=e()EFlfDHO)=FlY). 


Zur Ermittlung des Minimums wird die erste Ab- } 
leitung gebildet und gleich ‚Null‘ gesetzt. 


d 5 
a Ua+N LE = 


Fl) +0)f()=1 Flf() +0) = | 
dz 1 
a3 a 
an dz \da)zl 


Das Minimum wird also dann erreicht, wenn die 
Differentialquotienten an den Stellen x und z gleich 
sind, mit anderen Worten, wenn das betrachtete Flä- 
chenstück ‘in der Verteilungskurve von gleich großen 
Begrenzungsordinaten %’, und %’z begrenzt wird. 

Auf Grund dieser Tatsache wird nun ein graphi- 
sches Lösungsverfahren zur Ermittlung der günstig- 
sten Verteilungsbreite entwickelt. "Voraussetzung für 
die Durchführung des Verfahrens ist, daß bereits die 
‚ausgeglichene Verteilungskurve vorliegt. 

Zu dieser Kurve wird in bekannter Weise die Inte- 
gralkurve konstruiert. Damit der Rahmen der Zeich- 
nung nicht überschritten wird, wählt man als Einheit 
in der Ordinatenrichtung % oder 4, der Einheit der 
Abszissenachse. Hierzu wird noch eine dritte Kurve 
(Lösungskurve) in folgender Weise gezeichnet. Im 
Punkte %, wird die Senkrechte errichtet, die die Ver- 
teilungs- wie auch die Integralkurve schneidet (Bild 2). 
Durch den Schnittpunkt mit der Verteilungskurve Py, 
wird die Parallele zur &-Achse gezogen, welche die 
Verteilungskurve zum zweitenmal in Py, schneidet. 
Das Lot auf die &-Achse schneidet die Integralkurve 
in Pj,. Die Parallele zur x-Achse durch Py, trifft die 
Senkrechte über x, im Punkts Q. Der Natur der In- 
tegralkurve entsprechend ist der Abstand Q — Py, 
gleich dem Flächeninhalt der Verteilungskurve zwi- 
schen x, und &,.. Die Strecke Q — Py, wird nun von 


x 


7, abgetragen. Das Ver- 


N AURNAR: Bild 2, 


“ laufen läßt. Auf. diese Weise ergibt sich eine Reihe 


‚von Punkten P, P’, P’” usw., welche durch einen 


, stetigen Kurvenzug verbunden werden (Bild 3). Der 
auf diese Weise gewonnene Kurvenzug stellt die Be- 
ziehung F=f(x) dar. Definitionsgemäß sollte sich 
das neue Streuungsmaß auf die Fläche 4, beziehen. 
‚Eine Parallele zur x-Achse durch F/2 schneidet diese 
Kurve in R und die Integralkurve in 8. Das Lot von 
R auf die x-Achse schneidet die Verteilungskurve in 


Po. Die Parallele zur &-Achse durch De schneidet 


die Verteilungskurve zum zweitenmal in Py,. Die 
Differenz der beiden zugehörigen Argumente ©, —z, 
= v ist dıe gesuchte günstigste Verteilungsbreite. Das 
Lot von 8 auf die x-Achse gefällt, ergibt den Punkt W. 
Bekanntlicherweise ist dies der sogenannte Zentral- 
. oder wahrscheinlichste Wert, der die Fläche unter- 
halb der Verteilungskurve halbiert. 


\ 
| 


Bild 3, 


Die praktische Durchführung des oben dargestellten 
Verfahrens bietet keine besonderen Schwierigkeiten, 
der Zeitbedarf ist etwas größer als bei der rechnerischen 
Behandlung einer normalen Verteilung nach der üb- 
lichen Methode. Demgegenüber steht aber die Tat- 
sache, daß das dargelegte graphische Verfahren einer 
allgemeinen Anwendbarkeit zugänglich ist und nicht 
auf die sogenannten normalen Verteilungen beschränkt 
ist. Als einzige, wesentliche Schwierigkeit macht sich 
zuweilen die Vorbereitung der Verteilung für die Durch- 
führung des Verfahrens geltend. Voraussetzung ist 
nämlich eine gut.ausgeglichene Verteilungskurve, da- 
mit die gewonnenen Ergebnisse am Schluß nicht illu- 
sorischen Charakter tragen. Verfasser geht hierbei 
folgendermaßen vor: s 

Erfahrungsgemäß läßt sich bereits während der Er- 
hebung der Statistik mit großer Sicherheit eine Aus- 
sage über die Stabilität der Verteilung treffen. Zeigt 
sich hierbei im allgemeinen ein ruhiges und ausge- 
glichenes Verhalten, so werden die beim Abschluß er- 
haltenen Werte, als mit nur sehr geringen Zufälligkeiten 
behaftet, betrachtet und durch einen gleichmäßigen 
glatten Kurvenzug ‚‚nach dem Augenmaß‘ ausge- 
glichen. Zeigt sich jedoch bereits während der Er- 
hebung ein unruhiges Verhalten der Statistik, so daß 
sich am Schluß’eine recht zerklüftete Verteilungs- 
kurve ergibt, so empfiehlt sich eine sorgfältige Aus- 
gleichung nach dem sogenannten „Verfahren der glei- 


eben noch zuträglich ist. 


"Über . eine Zustandsgleichung mit Berück- 


tenden Durehschnitte“, Dieses Verfahren kann in 
tzt, indem man den Punkt P 


wiederholter Weise so oft auf die Verteilung ange- B 

wendet werden, als es für eine gute Ausgleichung not- 
wendig erscheint, bzw. als es für die hierdurch ver- 
ursachte Deformation der ursprünglichen Verteilung 


Haldensleben. h K. Wisser ot NER £ 


x 


sichtigung des kritischen Koeffizienten. 
“ Unter der Annahme, daß in einer homogenen im 
Gleichgewicht befindlichen Flüssigkeit der äußere 
Druck p[kg/m?] vermehrt um den Kohäsionsdruck 
gleich dem thermischen Abstoßungsdruck plus einer 
zusätzlichen (molekularen) Abstoßung ist, gelangt ei 

man zu folgender Zustandsgleichung: ee; 


aB® @ RT b 
+ 


+ 2 le een de Rn f 
Hierin bedeuten: ER B 
® in m?/kg das spezifische oder kg-Volumen, 
T die abs. Temperatur in Grad Kelvin, 
R die spezielle Gaskonstante in kgm/Grad, a 
a,b und n aus den kritischen Daten und der Gaskon- 
‚., „stanten R zu bestimmende Festwerte, a; 
‚pr [kg/m?], va [mP/kg]und Tx[°K] die kritischen Daten. 
Zur Bestimmung der 3 Festwerte a, b und n stehen 
3 Gleichungen zur Verfügung und zwar: 


1. die Zustandsgleichung, 


rn 
A 
eo 


2. die Bedingung: 22), BR 0, fürv — vg 
(horizontale Tangente im kritischen Punkt), 
3. die Bedin a a — 05% Bir — 
k gung: lauern, PR 
(kritischer Punkt = Wendepunkt). F 
Nach Einführung des kritischen Koeffizienten k = Seh. 
x Pk 2277 PkVk 


oder dessen Kehrwertes = 1/k = ergeben sich 


RTe 

die Festwerte von Gl. (1) wie folgt: 
a=pyvilds—), 

1 


1—2%' 


WE PRk(L— 27)? 


Bir 2(43—]) 


und n = 


Nach Einführung der dimensionslosen oder reduzierten 
Drucke, Volumina und Temperaturen: 


Pa=P|PR da = vlvor und T, = T/Tr 

erhält man die dimensionslose Form von Gl. (1): : 
Ti 1 (1— 22)? Ss 

Pa= ans (An—Ieit aa—ı)m Ab B 


in welcher außer den (dimensionslosen) Zustands- 
variablen ?z, v2, Tx nur der kritische Koeffizient k 
bzw. dessen Kehrwert & vorkommt. 

Die Kenngröße x bewegt sich etwa in den Grenzen 
von 0,2 (Essigsäure) bis 0,327 (Helium). Während sich 
aus bekannten Zustandsgleichungen nur ein einziger 
für alle Stoffe gültiger Wert von & (z. B. ®/, = 0,375 
bei van der Waals und */,, = 0,267 bei A. Wohl) ergibt, 
ist in der Gl. (2) der wirkliche, von Stoff zu Stoff ver- 
schiedene Wert des kritischen Koeffizienten verwendet, 

Für verschiedene Stoffe ergeben sich sonach ver- 
schiedene (spezielle) Zustandsgleichungen, wie z.B. 
für Kohlensäure (CO,) mit % = 0,276: 


Hy a 
BB tm nen. (3) 


128 N achrichten 


oder für Wasserstoff (H,) mit & = 0,306: 


Ta 4,465, 2,195 j: 
ER 22 + are Eh 


Für Stoffe mit x < 0,25 wird (4x — 1) negativ, und 
die Glieder 2 und 3 der rechten Seite von Gl. (2) ver- 
tauschen ihre Vorzeichen; so lautet z, B. Gl. (2) für 
Wasser mit x = 0,225: 


13,45 
use $ 


—— a 
Ug v2 vu ( ) 
oder für Ammoniak (NH,) mit, x = 0,243: 


Tg 35,7 . 38,82 
Bareln DET ec (6). 


Pz = 3,27 


Für große Volumina vz geht Gl. (2) in die dimensions- 
lose Gasgleichung: 


über, die mit pv= RT identisch ist.‘ 
Berücksichtigt man die Temperaturabhängigkeit des 


Kohäsions- und Abstoßungsgliedes in erster Näherung 
durch den Ansatz: ; 


ah N 
"2% T 4s—))e 


(1— 2x)? 


BE IETE 
FT aan (8) 


Pz 


so wird die Anfangsneigung der Dampfdruckkurve für 
den kritischen Punkt im 9, T,-Diagramm: ’ 


In Zahlentafel1 sind die nach Gl. (9) berechneten 
Werte von tg& den experimentell bestimmten (nach 
W. Schüle, Techn. Thermodynamik, II. Aufl. II. Bd. 
S. 80) für einige Stoffe gegenübergestellt. 


Zahlentafell. 


gem 2 17 27.28 "15,54 | 6,248| 7,2351 7,9 
tg& gemessen . .... . | 5,11 | 6,400| 6,910| 7,4. 


Die berechneten stimmen mit den gemessenen nur 


‚annähernd überein; doch ist der Gang beider Reihen 


der gleiche und der Vorteil ersichtlich, den die Ein- 
führung des kritischen Koeffizienten in die Zustands- 


. gleichung mit sich bringt. 


Wien. C. Kämmerer. 
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NACHRICHTEN 


Mathematiker-Tagung in Tübingen. 

In der Zeit vom 23. bis 27. September 1946 fand in 
Tübingen die erste Mathematiker-Tagung seit dem 
Ende des Krieges statt. Sie war von den dortigen 
Kollegen aufs sorgfältigste vorbereitet worden. Etwa 
140 Mathematiker aus allen Zonen Deutschlands und 
dem Ausland nahmen an ihr teil. Von den 44 Vor- 
trägen gehörten die folgenden fünf dem Gebiete der 
angewandten Mathematik an: 


O. Baier, Über die Hurwitzschen Bedingungen. 


L. Collatz, Lösung gewisser Differentialgleichungen 
mit dem harmonischen Analysator. 


G.Hamel, Aus der analytischen Mechanik. 

W. Quade, Auflösung linearer Gleichungen durch 
Matrizeniteration. 

G. Schulz, Das Summenproblem bei mehrdimen- 
sionalen arithmetischen Wahrscheinlichkeitsver- 
teilungen. 

In besonderen Sitzungen wurden ferner. organisato- 
rische, sowie Unterrichts- und Prüfungsfragen be- 
sprochen. 

G. Schulz (Aachen). ' 


2. -Angen Matti, Mech. Mur 
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